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SUMMARY 


Crystals growing from solutions with plane faces must have considerable 
supersaturation at some parts of the surface. A supersaturation number o 
is introduced to give a measure of the maximum supersaturation. For 
two-dimensional growth in very dilute solutions a general method for 
evaluating o is given; it is applied to the growth of regular polygons, to 
the growth of a step on a surface, and to growth on the edge of a plate. 
For three-dimensional growth only a rough estimate for the supersaturation 
occurring with regular polyhedra can be given; it is estimated that the 
maximum supersaturation on the surface of a cube or regular octahedron 
is between 0-25 and 0-4 of the concentration difference between the 


vertices and infinity. 


§1. [IyrRoDUCTION 


Iz is known that crystals can grow from solutions with plane faces. We 
consider here only those cases in which there is no particularly large 
surface mobility for solute atoms, and in which the solution is dilute, 
and the diffusion coefficient independent of concentration. Then at any 
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moment the flux of matter or the normal gradient of concentration has 
to be constant all over a face. This condition is the only one at the 
crystal surface that is required to specify the mathematical solution of 
the diffusion problem. This means that the concentration itself will not 
be uniform at a given face. All over a growing face, however, there must 
be at least a slight supersaturation; therefore a considerable super- 
saturation will occur on parts of the faces, especially near corners and 
vertices. This conclusion is in agreement with experimental observations 
on growth of square or rectangular crystals between plates from NaClO; 
solutions (Berg 1938, Bunn 1949, Humphreys-Owen 1949). It is 
presumably this excess of supersaturation at crystal corners which in 
certain circumstances causes dendritic rather than polyhedral growth : it 
is therefore very desirable to estimate how this supersaturation excess 
depends upon crystal geometry. 

The fact of growth in plane faces (plane means here only macroscopically 
plane, not plane on a molecular scale), in spite of a certain amount of 
supersaturation, is in accordance with the view (Burton, Cabrera and 
Frank 1951) that below a certain supersaturation crystals can grow only 
by the screw dislocation mechanism and not by a nucleation process at. 
the surface. 

This paper is concerned with theoretical investigation of the diffusion 
problem and the supersaturation conditions at the crystal surface. This 
would be a very difficult task for three-dimensional growth and for 
non-dilute solutions. The problem, however, is tractable for two- 
dimensional growth in a simple shape at very small concentrations. 
The latter condition ensures (together with the reasonable assumption 
that surface migration is unimportant) that the flux has to be constant. 
over any face, and that the growth rate is slow, so that the whole problem 
can be treated statically. This should give good results for the problem 
of the advancing step, because the characteristic length of this problem, 
the step-height a, does not exceed a few hundred angstréms, and is small 
compared with the ratio of diffusion-coefficient D and growth-rate v in 
a typical experiment. In the case, however, of regular polygons. 
(a=length of a side) growing in two dimensions it will be more difficult 
to satisfy the condition a<D/v, and our treatment is to be regarded as 
giving the limiting case for very dilute solutions. 

We cannot use any of the present observations on two-dimensional 
growth for comparison with the theory, because the conditions were not. 
sufficiently symmetric and the solutions were far from being dilute. 
However, it is hoped that future experiments at moderate concentrations 
and under well-defined conditions will supplement the theory in finding 
the general dependence of supersaturation conditions on concentration 
and growth rate. 

In §4 we give a few remarks on three-dimensional growth. They are 
meant only to give a very rough idea about the supersaturation conditions 
that are to be expected in three-dimensional growth of polyhedra. 


with the Growth of Crystals from Dilute Solution 3 


§2. GENERAL Metuop ror Two-Divensionan GRowTH 

The assumptions made in §1 permit the application of the methods 
of the theory of functions. The mathematical problem is to find the 
concentration ¢, satisfying the two-dimensional potential equation V2c—0, 
to which the diffusion equation reduces, subject to the flux conditions 
stated above. 

Two steps are involved in the solution of this problem. Firstly we 
consider the plane in which the two-dimensional crystal lies as a complex 
plane w=u-+iv, and we map the exterior of the crystal by a conformal 
transformation on the upper half of a complex plane z=x-+-iy. During 
this transformation the concentration c remains invariant, the normal 
gradient at the surface, however, is changed. Secondly we solve the 
potential problem in the x-y-plane by means of Fourier integrals. 

For a crystal with polygonal shape this can be done quite generally. 
The integrals involved can be evaluated only in simple cases. We are 
not, however, very much interested in the determination of the con- 
centration c(w, v) everywhere in the exterior of the crystal, but mainly 
in the supersaturation conditions at the crystal surface, especially in the 
largest supersaturation that occurs in a certain mode of growth. This 
can be characterized by a single number, the swpersaturation-number o, 
which we define as the difference between the highest and the lowest 
concentration at the crystal surface, divided by the integral of the normal 
gradient of the concentration, taken over the trace of the crystal surface 
in the w-plane. We shall evaluate o for several crystal shapes. 

According to the Schwarz—Christoffel transformation the interior of 
the polygon shown in fig. 1 in the w-plane is mapped on the upper half 
of the z-plane by a mapping function 


U2) Aen Mets SL. os ak Ga eRe) 
which satisfies 
Oe ple) =ple— 2,2... (22,4 BAe Se eres) 
or p=d(2)= [wo dee neared ate (3) 


Z1) «+++ 2m, being points on the real z-axis, are the images of the corners 
of the polygon. It is mainly the difficulty in finding the 2, for a given 
polygon that limits the present method either to very symmetrical cases 
or polygons with only few corners. ; 
The normal derivative (dc/dy),,-.-=G@(«) at the boundary in the z-plane 
can be expressed in terms of dc/dn and (x) according to the formula 
dc dcdn 
ol ee (4) 
dy onoy’ 
and is therefore known. The general solution of the potential equation 
V2c=0 with the right behaviour at infinity is given by 


ole, y= [_ g(keit* on | Mvdk. oie a ee 


BZ 
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From 
G(e)=— | gk) [leur dk se) ssa ein 6} 
g(k) is determined by 
il i : 
it=— srg | Gizye te de. 2 oe a 


The equations (3), (5) and (7) contain the solution of the problem, and 
practical applications are solely a question of evaluating these integrals. 


Fig. 1 


Notation for the angles of an arbitrary polygon. 


§3. EXAMPLES 
(i) Regular Polygons (fig. 2 (a)) 
By symmetry it is sufficient to consider the polygon ABCD with the 
interior angles %;=%,.—7($+1/n) if the regular polygon has » corners. 


The upper half of the z-plane (fig. 2 (b)) is mapped on the interior of 
ABCD by 


z 
w= | (L—#2)m—4) gp (8) 
, a, 
The factor 4 is to be determined from the fact that :——1 should 
correspond to w=—}a, z=1 should correspond to w= 3a (a being the 


length of a side of the regular polyhedron). Using the f 
(Groebner—Hofreiter 1950) i g the tormula 


1 
[e124 de BA, K=I(KIAI(eba, 2 2. (9) 


we obtain 


Lh) Pea 
a=nB (5,5 +3). eS ae fe ae) Rie ces. (10) 
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Fig. 2 


vr 
\ A ; 


(0) 
The region within ABCD of fig. 2 (a) is mapped upon 
the region y>0 of fig. 2 (0). 
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Let the concentration gradient at the boundary of the regular polyhedron 
be dc/@n=q. Then the boundary conditions in the z-plane are given by 


0  |a|[>1, } (11) 
G2 a aes ape’ "otha 


Due to the symmetry of G(x) we may write the general solution for 
C(x, y) as 


o(u, y)=2 | g(k) cos kee e-¥Y dl, on i tee RGL ot 
0 
Where g(k)=— 5, |” G() cos ke dé. ee aaa) 


It is a mathematical convenience to introduce the complex function ¢(z), 
whose real part is the physically significant concentration c(z, y)=Zc(z). 


Then czy=— = [" [SBE tect) ae ak. Cy ea 
/ 7 0 0 


In order to obtain known functions we consider de¢(z)/dz, interchange the 
order of integration, and use the formula 


= is ike —— wie 1 
it cos ke et dk=— a, . (15) 
which holds for y>0. 
This gives 
dc(z) DRT z 
os ——— o) tne gee ea 
ae (16) 
and (apart from an unessential additive constant) 
| r 2 = 
e(e)= = | G(é) In (@—é) dé. “| > Se 
0 
Using eqn. (11), eqn. (16) becomes 
dc(z)___ 2qu (1 2)(1/n—1]2) _* 
=— Fl aa apd =. . (18) 
This integral can be evaluated in terms of the hypergeometric function 
pe Ic) : b-1 e—b—1 —a 
F(a, b, ¢, 2)= Teresa |! (1—#)e-?-(1—tz)-4 dt... (19) 
The result is 
de(z) qu ,/l 1 il lel 


By analytic continuation this holds not only for y>0, but also for y=0. 
We did not succeed in expressing ¢(z) in terms of well-investigated 
functions. It is, however, possible to express the minimum concentration 
c(0) at w=0 and the maximum concentration ¢(1) at w=ta in terms of 
the logarithmic derivative of the factorial 


— dinii(z) din I(z+1) 


P(z) Ae er ea ake (21) 


with the Growth of Crystals from Dilute Solution 7 


by using the formula (Groebner—Hofreiter 1950) 


1 
| x-1(1—#)*1 In (1—a) de = B(x, )[P(A—1)—P(e-A—D)]. 


0 
Bei) -O) 


This gives 
_ be Ui. Sah 1 1 1 
The saturation number is 


= S530 [e(-34)-r(-3)) « 


An approximately valid formula for large n is obtained from the leading 
term, 


and 


o=7/4n?, ee eae eee ee 26) 

of the Taylor expansion in powers of 1/n. It gives values for o (in table 1 

given as o,) which are too large, the correction term being of the order 1/n’. 

Another upper bound for o was suggested by F. C. Frank. If the 

concentration field were that of a line-sink of unit strength (so having 

circular symmetry despite the polygonal form of the crystal boundary), 
giving 


(Peay PZ AM NON Dae EG aes en Yelerg AAD) 
the supersaturation number would be 
o=—(1/27)In cos(aj/m). «|. . » - = (28) 


The actual value for o (which we give in table 1 as o,) must be between 0 
and the value given by eqn. (28). For large n eqn. (28) reduces to 
eqn. (26), the next term being of the order 1/n’. 

We refer to the approximation that gives eqn. (28) as the ‘ circle- 
approximation ’; corresponding values of o appear in table 1 as o;. The 
corresponding ‘sphere-approximation’ may be applied to regular 
polyhedra and will be dealt with in §4. ; 


Table I 


0-111 | 0-0570 | 0-0349 | 0-0171 | 0:0101 | 0-00674 | 0-00177 


oo 0-1103 | 0-0551 | 0-0229 | 0-0126 | 0:00798 | 0-00198 
0-196 | 0-0822 | 0-0492 | 0-0218 | 0:0123 | 0-00785 | 0-00196 


As it is difficult to give the concentration along one of the faces of the 
polygons by evaluation of eqns. (8) or (17), we shall now derive a simple 
interpolation formula for the concentration along the faces. We denote 
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this concentration as a function of u by C(u), and we express O(w) by an 
interpolation formula giving the correct value for o and for dC (u)/du at 
u=1ta (which we denote by C’(3a)). Such an interpolation formula is 
C(u)—C(0) = {20 (3a) — 20(0) — a0 (4a) }(2uJa) 
+ {4a0"(3a)—C(3a)+C(0) }(2u/a)*. (29) 

In order to evaluate C’(}a) we apply to eqn. (20) the formula 
I'(c)I'(c—a—b) 
I'(c—a)I'(c—b) 
I(c)I'(a+b—c) 
__»\c—a—b 

ee Seelayit) 
and proceed then to the limit w=}a(z=1) in the equation for C’(u). 
This gives after some calculation 


Dhaes a ie ea 
r(5+; r(5—; 
vq ee eee es! sie een 


T I(1) cos (z/n) 


Our interpolation formula (31) becomes now 


ee {20n— restate) + {rent eT Ge). 
(32) 


This method can of course be extended to include higher derivatives, for 
instance d°C/du* at the centre of an edge. 


Pelee C1 ae F(a, b; a+b—c+1 ; 1—z) 


F(c—a, c—b ; c—a—b+1; 1—z) . (30) 


(ii) Growth of Steps on Crystal Surfaces 

Crystals showing growth-spirals grow in steps of a step-height equal to 
the Burgers vector of the screw-dislocation which forms the axes of the 
spiral. One does not know anything about the actual shape of the step. 
We assume here that it is a plane, but we do not make any assumptions 
about the angle y (see fig. 3). We can treat this problem along the lines 
indicated in §2. We introduce a complex variable w=u-+iv as shown 
in fig. 3. The mapping function of the (degenerate) polygon formed by 
the trace of the crystal surface on the z-plane is 


Zz 1-3 yln : 
w=p| (=) ine eet ene eae 
0 
From the condition that the distance between A and B is equal to a, 
is determined by 
=, [> i( dhe eee 
a= (5) dsm eect tenitea) 


The function (dc/dy),_»—G(«) is, if the normal gradient of the concen- 
tration along the step is denoted by 4, 


0 (ze eel \ 


G(x)= 
im | tee |e | <1. be 
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By calculations similar to those in the first example one finds for the 
concentration 


ce, y= Bey =" @ fe me—e (ES) "ae . . (36) 


I—¢ 
The difference 4c between the concentrations at z—=—1 and z—1 is given 
by 
Baie Hae \ Te I 
Go a iS (=) In (=) dé=qa (; — cot ”) pee eae CG 7.) 
and the saturation number is 
_Ac_ i taney ae l eee 
ram wo a3) an? gag? 1 °° gee ASS) 


o is plotted as a function of y in fig. 4. A remarkable feature is that o is 
approximately a linear function of y for small y but that it rises sharply 
if the angle becomes an obtuse one. 

We should mention that the problem solved in this example also 
provides the solution for a crystal growing in the shape shown in fig. 5. 


Fig. 3 


w=utiv 
ep 


| 
| 


Coordinates for a stepped crystal surface. 


§4. THREE-DIMENSIONAL GROWTH OF REGULAR POLYHEDRA 


Quantitative observations during growth of crystals suitable for 
checking the theory are most likely te be made for two-dimensional 
growth. For natural growth of crystals, however, the growth in three 
dimensions is most important. The simplest cases of growth in polyhedral 
shape are growth in the form of a cube or of a regular octahedron. These 
cases cannot be treated exactly even under the idealized conditions which 
we introduced into the earlier parts of this paper. These idealizations 
allow one to replace the diffusion problem by a problem. of heat con- 
duction ; we have to find the temperature outside a solid body with zero 
thermal conductivity and with a uniform distribution of heat sinks on 
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Fig. 4 

Oo 

fs) 

1-0 

0:5 

Tt/4 It/2 3% Y 
Dependence of the supersaturation range at a step upon its steepness, y. 
Fig. 5 


g 
i 
| 
| 
| 
I 


Another crystal form to which the step solution applies. 
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its surface. In this heat-conduction problem we can ask for the solution 
of this problem not only for octahedron, cube and tetrahedron (these 
are the only regular solids that are possible growth shapes of crystals), 
but also for icosahedron, dodecahedron and for the dihedra. These 
dihedra (Klein 1924) are improper regular solids with zero volume ; they 
consist of two regular polygons as the only faces. 

We shall apply the sphere-approximation to the proper regular poly- 
hedra, and we may expect from the results in §3 that this will give an 
overestimate for o of about 10 to 20% for the icosahedron and the 
_ dodecahedron. On the other hand we can solve the heat-conduction 
problem accurately in a simple way for one of the dihedra, namely the 
limiting case of a large number of corners, the circular disk. So we get 
for the sequence of regular solids two numerical values for the saturation 
number comparatively accurately, and we may interpolate for the other 
regular solids. (We continue to give the discussion in terms of 
concentration rather than in terms of temperature.) 

We denote the radius of the circumscribed sphere of the regular 
polyhedra by R,, the radius of the inscribed sphere by &,, and the length 
of the edges by d. In order to get a dimensionless quantity for the 
supersaturation number o we adopt a definition slightly different from 
the one used in two dimensions. We define o/R, as the difference between 
the maximum and the minimum concentration divided by the integral 
of the normal gradient of concentration, taken over the whole surface. 

The concentration field of a unit point-sink is given by 


Ca Aa ee a) se. Ae on (39) 


In the sphere-approximation we assume this to be the concentration 
also in the presence of the solid. The supersaturation number is then 
given by 

rer (hida\ ice nina ee ae nel ae (40) 


In table 2 we give a list for R,, Ry (Buerklen—Ringleb 1936) and for o, 
calculated from eqn. (40). 

For the circular disk the sphere-approximation would give an infinite 
supersaturation number. We can, however, solve this problem exactly. 
If p, z, ¢ are cylindrical polar coordinates, p=R,, z=90 correspond to the 
circumference of the disk, and ¢ is equal to | dc/dz | on the disk, then the 
concentration is given by 


c(p,2)=—aR, J, Foe) TAR,)er* dr). mide boats (AL 


J. is the Bessel function of order v. The minimum concentration at the 
centre of the disk is 


(0, 0)=—4R, | TORN =the oe) 
0 
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the maximum concentration at the circumference is 
o(R,, 0)=—qR; | Jo(AR,)F,(AB,) dAJA=—2QR,/7. . . (43) 
0 


The supersaturation number is 


1 1 

c— on ee = (005795 Bo eet eee (44) 
This is less than the value of o in table 2 for the tetrahedron, cube and 
octahedron. Comparing with n=2 in two dimensions one would expect 
the proper solids to have smaller supersaturation numbers than an 
improper one. This shows that we have to reduce these values of o to 
a considerable extent. We think that it might be a good estimate to 
assume that o is between 0-025 or 0-03 for cube and octahedron, and 
between 0-03 and 0-04 for the tetrahedron. 


Table 2 


R,/d R,/R, 


[(7+3/5)/6] | 1,258 
[(50+224/5)/5] | 1,258 
6 1,732 
1,732 
2,000 


Icosahedron 
Dodecahedron 
Octahedron 
Cube 
Tetrahedron 


ay 
4 
at 
4 
i 
2 
1 
2 
Al 
4 


In three dimensions it is possible to relate the concentration difference 
A,c between infinity and the surface of the crystal to the supersaturation 
Ac at the surface. If the concentration field were really that given by 
eqn. (39), and if 4,¢ were taken at the vertices of the solid, this relation 
would be 


Ac=4ra Ac. 5. Para cee ooo eee) 
In the case of the disk the corresponding relation is 

c= 057110, a nT: 
if 4,¢ is taken between the edge of the disk and infinity, and 

c= 0-364 Arh mee eee SS en aaa yy 


if 4,c is taken between the centre of the disk and infinity. 

For the disk the numerical factor 470 in eqn. (45) would be 1-23 times 
the correct factor given by eqn. (46), so that both give the same order 
of magnitude. A rough estimate for a cube and a regular octahedron 


would be that under symmetrical conditions dc is b 
é etween 25° 
40%, of Aic, ea 
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[Received October 15, 1952] 


ABSTRACT 


Further studies of the helium film have been made using the optical 
method of Burge and Jackson. Evidence has been obtained that there 
is a film of helium I at 2-3° Kx about 10 atoms thick, but accurate values 
of this thickness have not yet been obtained. It has been shown that 
the thickness of the creeping film of helium II does not change when 
the volume of liquid flowing per second in the film decreases considerably 
and that the thickness is equal to that of the stationary film under the 
same conditions. 

Bulk liquid is seen to form when the helium II film flows on to a 
surface (i) which is below the highest liquid level and (ii) which has a 
perimeter smaller than that which controls the rate of creep. The bulk 
liquid returns to the film if the latter then flows on to a surface with a. 
perimeter at least as great as the controlling perimeter. 


§1. InTRODUCTION 


In a previous paper (Jackson and Henshaw 1950)T a description was 
given of simultaneous measurements of the transfer rate of liquid 
helium IT from a stainless steel beaker and of the thickness of the moving 
film on the outside of the beaker by the optical method of Burge and 
Jackson (1951). These observations have been continued for the 
investigation of certain specific problems. A preliminary account of 
part of the work was given at the Oxford Conference on Low Temperature 
Physics, 1951 (Jackson and Henshaw (1951)). 


§2. OBSERVATIONS oF Hetium I Fim 


THE work of Kistemaker (1947) and Long and Meyer (1949) has shown 
that one may expect to observe a film 15-30 atomic layers thick on any 
solid substance in contact with liquid helium I at temperatures 
immediately above the \-point. Such a film was not observed by Burge 
and Jackson but their technique was, however, as follows. Readings of 


* On leave of absence from the National Research Council of Canada. Now 
at Atomic Energy of Canada Limited, Chalk River, Ontario. 
+ Communicated by the Authors. 
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the position of the nicol prism were made at some temperatures below 
the A-point. The temperature was then allowed to rise above the A-point 
and was maintained constant at, say 2-3° x. Readings of the position 
of the nicol prism were again made, after which the radiation from a 
small electric lamp held close to the cryostat was allowed to fall on the 
mirror for a few minutes. After switching off the lamp the position of 
the nicol prism was read once more. The averages of the readings of the 
second and third sets were always found to be the same within the 
accuracy of the observations. This could mean (i) that no film was 
present at the temperature 2-3° K or (ii) that a film was initially present 
at 2:3°K and was evaporated by absorption of radiation from the lamp 
but was reformed completely either from the bulk liquid or from the 
vapour phase, after the lamp had been switched off, during the period 
of 1-2 minutes required to take the first few observations of the third 
set. The first alternative was assumed by Burge and Jackson and the 
difference between the first and second sets of nicol readings was taken 
as a measure of the total thickness of film present at the temperature of 
the first set of observations. 

The second alternative has now been investigated and found to be 
the correct one. For this purpose a small electric heater was attached 
to the top of the mirror used by Burge and Jackson, the current leads 
being taken down into the liquid helium and then up the support. Thus 
no heat was conducted down the leads directly to the mirror and the 
helium film on it. Observations were then made of the nicol position 
(i) below the A-point (no heat input), (ii) at 2-3° K (no heat input) and 
(iti) at 2-3°K (heat on). For a dissipation of 50 milliwatts or more the 
nicol setting was seen to change from (ii) and take on a new constant 
value. On switching off the current in the heater the setting returned to 
the value (ii) within about one minute. The nicol shift leet 
settings (ii) and (iii) corresponded to the presence of a film at 2-3° K 
of the order of magnitude of 10 atomic layers* and so all film thicknesses 
given by Burge and Jackson (1951) and Henshaw and Jackson (1950)f 
have to be increased by this amount. Owing to circumstances accidentally 
present, these observations could not be made with the accuracy usually 
obtained. We shall therefore not give any table of corrected values here 
but shall postpone publication until the results of definitive experiments 

i ress are available. 
Be Caisalitions of the velocity of the helium II film as given by J. ackson 
and Henshaw (1950, 1951) involve not only the measured film thickness 
but also the measured creep rate and the geometric internal perimeter of 
the beaker. The recent observations of Chandrasekhar and Mendelssohn 


* Taking the value 4-0 x 107° cm (Keesom and Schweers) for the thickness 


ic layer, as used by Kistemaker (1947). 
Se oa tiiclnssd measurements quoted by Henshaw and J ackson at the 
Conference on Low Temperature Physics at the National Bureau of Standards, 


Washington, 1951. 


16 L. C. Jackson and D. G. Henshaw on Further Observations 


(1952) with highly polished stainless steel appear to show that the higher 
creep rates observed on less highly polished metals under clean conditions 
are a mixture of true film flow at a universal rate depending only on the 
temperature and a flow in cracks and other channels. Until the present 
controversy as to the interpretation of the observed creep rates 1s finally 
settled it is perhaps inopportune to publish values of the calculated film 
velocities. 


§3. THICKNESS OF THE TRANSFER Fito AT SUBCRITICAL VELOCITIES 

In the experiments of Jackson and Henshaw (1950) in which the thick- 
ness of the helium film was measured during transfer into or out of a 
cylindrical stainless steel beaker, the observations were made on the film 
on the outside of the beaker. The liquid helium in the film at the point 


Fig. | 


of observation was necessarily moving with a subcritical velocity. It was, 
however, not known whether variation of the velocity within the 
subcritical range had any influence on the thickness of the moving film. 
To decide this question and to obtain visual evidence of the production 
of bulk liquid as described in the next section, the beaker shown in 
fig. 1 (a) was constructed. It was again of stainless steel and consisted 
of two cylindrical portions, the upper of internal diameter 0-7 em and 
external diameter 0-8 cm, the lower of 0-3 em and 0-4 em respectively. 
It was provided with a glass capillary tube to observe the level of the 
liquid in the beaker and with a radiation shield as shown in fig. 1 of the 
previous paper. A strip of plane mirror, 2 mm wide, was ground and 
polished on the outside of each portion of the beaker and barium stearate 


films were deposited as usual so that there was a 1-3 boundary at A and 
a 3-5 boundary at B. 
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Observations were then made at the 1-3 boundary at various 
temperatures below the )-point after the beaker had been filled and 
withdrawn almost entirely from the liquid helium. The thickness of the 
film at A was found to be the same within the accuracy of the observations 
when the inner liquid level was in the upper wider part of the beaker or 
in the lower narrower part and was equal to that of the stationary film 
present when the transfer had ceased. Thus since the volume of liquid 
transferred per. second decreased by a factor of more than 2 when the 
inner level entered the narrow part, the velocity of the helium in the 
film must have decreased by the same factor. For subcritical velocities 
of flow the film thickness therefore depends only on the temperature, the 
height of the point of observation above the main liquid surface and 
(presumably) on the nature of the substrate. The velocity of the flowing 
helium adapts itself to the amount of helium transferred per second. 


§4. PropucTIon or BuLK Liqurp FRoM THE FILM AND VICE VERSA 


When observations were made at point B on the lower part of the beaker 
the results were as follows. When the inner liquid level was in the wider 
part of the beaker the appearance at B was that of a waterfall or a stream 
of liquid with large local variations of thickness. When the inner level 
was in the narrower part the usual film was observed with a thickness 
corresponding to the height h,. Bulk liquid was therefore formed when 
the film carrying a volume of liquid per second proportional to the internal 
diameter of the upper part of the beaker flowed on to the narrower part 
of the beaker. Since however the thickness of the film at A was that 
for the height h, and not the much smaller height h, (the distance to the 
point at which bulk liquid first appeared) the film must have been 
continuous over the whole beaker and the bulk liquid represented the 
excess liquid which could not be carried by the film on the surface of 
reduced perimeter. 

Mendelssohn and Chandrasekhar (Mendelssohn (1951), Chandrasekhar 
(1952)) have studied the formation of bulk liquid from the helium II film 
with a glass beaker of varying internal diameter and provided with 
external ‘skirts’. They interpreted their observations as showing that 
bulk liquid is formed from the film at points below the level of the liquid 
in the beaker when the film flows on to a surface with a perimeter smaller 
than that governing the rate of transfer and that film is reformed if the 
total liquid flowing then encounters a surface of perimeter larger than 
(or equal to) that originally governing the rate of transfer. 

The above observations by the optical method give a complete visual 
confirmation of the correctness of the first part of their conclusions. To 
test the second part a further beaker of the shape shown in fig. 1 (6) 
was constructed.* The observations on the,upper and middle portions 
were exactly as recorded for the beaker of fig. 1 (a) as was to be expected. 
Observations on the lower portion of the new beaker showed the presence 
eee SS 

+ We are indebted to Mr. A. C. Ham for these observations. 
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of a film of a thickness corresponding to the height of the lower boundary 
above the outer liquid level independent of the position of the inner 
liquid level. The helium, which appeared as bulk liquid on the narrow 
part of the beaker, had therefore returned to the film form on the wider 
lower part of the beaker. 

A further simple observation* is of interest. A beaker was filled with 
helium II and partly withdrawn from the liquid. An accidental dis-~ 
turbance to the support caused some of the contents to splash over. 
This liquid was seen to run down the mirror in the form of separate drops 
(diameter 1-4 mm). In spite of the fact that the beaker was covered 
with a thick film of helium II the drops showed no tendency to spread. 
This may seem strange at first sight but it is what must be expected 
when the surface is covered by the creeping film of helium II which can 
carry only a certain maximum flow at any temperature. Liquid in excess 
of this flow cannot enter the film and remains as bulk liquid in the form 
of separate drops. When the observations were repeated at a temperature 
above the A-point the results were quite different. The liquid splashing 
over was seen to drain down and spread sideways in the manner to be 
expected with any normal viscous liquid which wets the substrate. 
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ABSTRACT 


The effect of the polarization field on the dielectric constant of a 
medium can be taken into account in two alternative ways, which are 
equivalent. It may be regarded as enhancing the actual field that 
produces the polarization, or as determining the resonance frequencies 
of the medium, which will be different from the frequencies of the 
individual oscillators that constitute the medium. The dispersions of 
the alkali halides are discussed from the latter point of view. Though the 
polarization factors defining the interactions between the infra-red oscil- 
lators have the Lorentz value 47/3, those defining the interactions of 
the electronic oscillators are found to be zero. The electronic overlap 
between the neighbouring ions in their equilibrium positions seems to 
be just the one that corresponds to negligible interactions between the 
electronic oscillators. 


$1. INTRODUCTION 


THE part played by the polarization field in dielectric phenomena has 
been discussed in detail by Lorentz (1908), and following him by several 
authors. (For a good account of the subject see Debye 1928, Fréhlich 
1949, or Rosenfeld 1951.) An obvious method of considering the influence 
of the polarization field on the dielectric constant is to regard the effective 
field that produces the polarization as including the polarization field in 
addition to ‘ the field in the medium ’, to which the polarization field will 
naturally be proportional. The result will be an enhanced susceptibility 
and a correspondingly enhanced dielectric constant. This method of 
calculating the effect of the polarization field is the one usually adopted. 
The polarization field is also known (Krishnan and Roy 1951) to play 
an important part in determining the observed resonance frequencies 
of the medium. Hence its effect on the dielectric constant may also be 
regarded as exercised indirectly through its influence on the resonance 
frequencies. The dispersion of the dielectric constant of the alkali 
halide crystals is discussed in the present paper from the latter point of 
view. The alkali halides are particularly suitable for this purpose. 
cs ES a 
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The point of view adopted here throws also some light on the relation 
between the different types of dispersion formulae, and in particular on 
the relation between the characteristic frequencies that appear in these 
formulae. 


§2. Tae Errecr or THE POLARIZATION FIELD ON THE PRINCIPAL 
Larrick OSCILLATION OF AN ALKALI HALIDE CRYSTAL 


In a recent paper (Krishnan and Roy 1951) the frequency of the 
principal mode of oscillation of an alkali halide crystal, in which the 
lattice of the alkali ions oscillates with respect to the lattice of the halide 
ions, was calculated on the basis of the simple Born model. It was found 
that the electric polarization of the crystal due to the separation of the ~ 
charges accompanying the oscillation, plays an important part in deter- 
mining the frequency. Separating the two lattices by a small distance r 
one can obviously express the potential energy of the crystal thus deformed 
as a power series in 7, in which owing to the centre of symmetry of the 
lattice points of the undeformed crystal, terms involving odd powers of 
r will naturally be absent. We shall confine attention at present to the 
7 term, which is the first non-vanishing, and also the predominant, term 
involving 7, and which determines the frequency. We shall denote this 
term in the potential energy per pair of ions by 


W=arr=thyw fr, ary ik (1) 
where yz is the reduced mass of the pair given by 1/u=1/m,+1/m, and 
m, and m, are the masses of the two ions, and w; is the angular frequency 
of the principal mode of oscillation. 
The essential part of the problem is to evaluate the contributions to 
a from the different types of interactions. This has been done in the 
paper referred to, and we shall merely quote here the relevant results. 
It will be seen immediately, in view of the centre of symmetry of all the 
lattice points, that the electrostatic interactions contribute nothing to a. 
As regards the repulsion interactions, they may be taken to be of the 
form Ae~*’, where R is the distance between the interacting ions, and 
they are of short range, since p is found to be much shorter than the 
equilibrium distance d between the neighbouring ions. Hence these inter- 
actions may be regarded as confined to the immediate neighbours only. 
Their contribution to a is found to be equal to 
3 ae?(5 — 2) 
ONG dt Gas met BAY, slut tay oe A ee 


where N is the number of ion-pairs per unit volume, B is the compress- 
ibility of the crystal, « is the Madelung constant, d is the equilibrium 
distance between neighbouring ions, 8=d/p and e is the electronic charge. 
The interactions of the van der Waals type involving the dipole moments 


induced in the ions contribute relatively little to a, and they are neglected 
throughout. 
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Had a, been the whole of the contribution to a, the frequency of the 
principal lattice oscillation of the crystal would be given by 


a= uQ 2. . ° ° . . . . . : (3) 
But actually the contribution to a from the polarization field is consider- 


able. For a given relative displacement r between the two lattices the 
polarization per unit volume is given by 


do CraR Ir Mea aot wi «a sce cus fot (4) 


In the paper referred to it was shown (1) that owing to the electronic 
overlap between the neighbouring ions, this polarization does not induce 
any electronic polarization in the ions, and (2) that the polarization field 
acting on an ion and tending to displace it as a whole, has just the Lorentz 
value, namely (47/3)P, and hence (3) that the corresponding contribution 
to a is given by 

Ao=—irNe*. RUN RS ets ie, rave Pe Lied Cs) 


This term is negative since the polarization field is in a direction that 
tends to increase the relative displacement r between the two lattices. 

In NaCl, for example, a, is found to be numerically nearly half of a, 
and since they are of opposite signs, a, is of almost the same magnitude 
as a, the net coefficient of the r? term, namely 


A=0,+4. SPT ee Debian Sik uses os CO) 

For most of the alkali halide crystals the frequency w, calculated in 

this manner by taking the polarization field into account, agrees well 
with the observed reststrahlen frequency (Krishnan and Roy 1951). 


§3. THE MAGNITUDE OF THE POLARIZATION FIELD 


As is well known the internal field F in a dielectric may differ consider- 
ably from the field Z in the medium, the difference between them being 


‘the polarization field . 
F=E+pyH, a 4 ° ° . e 5 . (7) 


in which y is the susceptibility of the medium, ie., the polarization per 
unit volume per unit field in the medium, and’p is the polarization factor, 
which is rather difficult to evaluate except under certain simple conditions. 
When the elementary dipoles in the medium, which togother make up the 
polarization of the medium, are point-dipoles, and when they are arranged. 
in a cubic lattice, or are distributed at random, p, as Lorentz has shown, 
will be equal to 47/3 (see also Van Vleck 1937). Pb. 

But actually, as was shown in the paper referred to, the dipoles concerned 
in refraction phenomena, for which Lorentz’s theory was originally 
intended to apply, are by no means point-dipoles ; they are due primarily 
to the displacement of the electrons in an atom with respect to its 
positively charged nucleus, and the separation of charges will be of 
atomic dimensions, and comparable with the distances between neigh- 
bouring dipoles. Thus the Lorentz polarization factor p=47/3 may fail 
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just where it was originally intended to apply, namely to electronic 
polarization as is involved in refraction phenomena in the visible and in 
the ultra-violet regions of the spectrum. But the polarization produced. 
in ionic crystals like the alkali halides by a small relative displacement of 
the two lattices, with which we are concerned here, is in a different 
category altogether, since we can now regard the polarization as due to 
dipoles of small dimensions of the order of 7 located at the lattice points. 
This can be readily seen from the following considerations. Electrostatic- 
ally the effect of deforming the original lattice, which we shall designate 
by A, into B will be equivalent to superposing on A the following two 
lattices: the lattice A’ obtained from A by changing the signs of all the 
ions in it, and the lattice B. The superposition of A’ and B will be equi- 
valent to locating a small dipole of moment er at each of the positive, 
or each of the negative lattice points, which again will be equivalent to 
locating a dipole e7/2 at every lattice point, whether positive or negative. 
Since r can be made as small as we desire, these dipoles will be practically 
point-dipoles, and the corresponding polarization factor p will therefore 
be expected to approximate closely to the Lorentz value 47/3, and this 
is just the value which we have used in the caiculation of w; and which, 
as we mentioned, is confirmed experimentally. 


$4. THe RESONANCE FREQUENCY OF AN ASSEMBLAGE OF 
HARMONIC OSCILLATORS 

Coming back to the calculation of the energy ar” per pair of ions from 
which the resonance frequency w,; of the crystal is evaluated using 
relation (1), it may appear at first sight that when the electrostatic 
interactions between all the ions of the positive and the negative lattices 
displaced relatively to each other by r are taken into account, the effect 
of the polarization field also, which is a direct result of this displacement, 
gets included in it automatically, and it should not therefore be necessary 
to invoke separately the effect of the polarization field. But actually 
the necessity arises from the following circumstance. It was shown just 
now that separating the two lattices by a small distance r is equivalent 
to developing small dipoles at the lattice points. The oscillating crystal 
will be equivalent to harmonic oscillators, all of the same moment, located 
one at each of the lattice points of the crystal. The frequency Q; which 
we would have obtained from a, with the help of relation (3), i.e. the 
frequency which we would have obtained by ignoring the polarization 
field, would be the frequency of these individual oscillators, whereas the 
frequency w; that we are trying to calculate is the resonance frequency 
of the assemblage of these oscillators, which will be different from Q, 
and the difference is determined by the polarization field. (That 2; is 
itself dependent on the density of the medium does not affect the argument.) 

Anticipating the results concerning the energies associated with the 
aco of the medium to be given in § 6, we may mention immediately 


9 
Nar? =3Nuw 2°=4PH, 
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where H is the field in the medium, whereas the energy that is obtained 
from the Coulomb and the repulsion interactions is equal to PF'/2, where 
Ff is the internal field, which includes in addition to E the polarization 
field pP. Hence one has to subtract from Na,r? the work done by the 
polarization field, namely pP?/2 (=27P2/3 in the case of the alkali 
halides), before it can be equated to 4Nyuw ?r?. 


§5. THe DIELEcTRIc ConsTaANT OF THE CRYSTAL AND THE 
DRUDE AND THE LoRENTzZ FoRMULAE 


Let us for the present ignore the electronic polarization, i.e. the 
polarization due to the displacements of the electrons in the ions with 
reference to their respective nuclei, and regard the medium as consisting 
of ions that can only be bodily displaced. For any given relative dis- 
placement r of the positive and the negative lattices the polarization 
per unit volume will be given by (4), namely P=Ner. This polarization 
may also be regarded as produced by a field H in the medium, such that 


eas 20a Ory eM Be ere) sat Ree. (8) 
The corresponding internal field F will be given by 
+ 
Fe=He+ 5 Netr=2ayr—pQ?r, se eae a) 


where, as we have seen, ; is the frequency which we would have obtained 
for the principal oscillation of the crystal if we had neglected the effect 
of the polarization field, as distinguished from the actual resonance 
frequency ,. 

It will be seen from (8) that the susceptibility x) of the medium for 
static electric fields will be given by 


Gs 
ig Ee ory . . . . . . . (10) 
where ‘ 
Be OM ee es nto LI) 
ae 


and ¢, is the corresponding dielectric constant. Equation (10) can be 
readily identified with the well-known expression of Born’s (see Handbuch 
der Physik). The two frequencies 2; and w, are connected by the 
relation 


Ge 
Gp ne ue etic (2) 
or 
2cot? 13) 
220 7 3 a . ° ° . . . . ( 
Similarly 
ilies 
pe a? ae (14) 
oa 
Ag=—a a 
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Hence eqn. (10) can also be written in the form 
é—l €,/8 
Ey +2 = QP? : 
Similarly the dielectric constant at any frequency can be expressed 
either in the form 


(15) 


e,—l= EP 


or, in view of (12), in the form 

€,—1 C,/3 

eee -e Q7—w’ 
which are indeed identical relations (Livens 1912). 

Just as the equivalence of (10) and (15) is based on relation (13), the 

equivalence of (16) and (17) corresponds to a relation which is analogous 
to (13), and which is a corollary to it, when we are concerned, as we are 
now, with only one frequency that influences the polarization, the 
relation being 


Bee ee 


Eg +2 

ge 
It can be readily seen that (16) and (17) reduce to (10) and (15) respectively 
when w— 0. } 

In (16) resonance occurs obviously when w=w, and in (17) when the 
right hand side tends to unity, i.e. when 2,7—w?=C,/3 which in view of 
(12) corresponds again to w=w, as it should. 

One may thus take into account the effect of the polarization field on 
the dielectric constant explicitly by taking the field that is effective in 
producing the polarization as (1+-47y,/3) times H or (e,+2)H/3 where 
E is the field in the medium, in which case the formula obtained is of 
the Lorentz type, and the characteristic frequency appearing in the 
expression is the frequency Q;; or alternatively, one may regard the effect 
of the polarization field as confined to changing Q; to the actual resonance 
frequency of the crystal w,; such that the relation (18) is satisfied, in which 
case the formula is of the Drude type. 


Q2—w? =(w2—«") 


(18) 


§6. Tue Evecrric Frerps Associatep witH THE SEPARATION 
OF THE LATTICES 


Now the potential energy Nar® per unit volume, due to separating the 


two lattices by a small distance r can also be regarded as due to the 


incidence of a ‘ field H in the medium’ such that eqn. (8) is satisfied, or 
such that 


=P. nce 9 Le 
Similarly ne) 


Na,r=tFP, ea Ae Sw OS greet) 


and the Resonance Frequencies of the Alkali Halide Crystals 25 


and the energy per unit volume associated with the occurrence of the 
polarization field will be given by 


Naor =—=tpP?, eae ae te cl) 
where pP is the polarization field. In the alkali halides, as we have 
seen, p=47/3. 

As is well known, the field in the medium £ at any given lattice point O 
includes not only the field due to all the surrounding dipoles, but also a 
certain averaged field at O due to the dipole at O also, i.e., the averaged 
self-field of the dipole at O also, whereas the inner field F' does not include 
the latter (Frohlich 1949). Though we refer to the latter field for 
convenience as the self-field, it depends very much on the polarization 
of the surroundings, which will be realized immediately when we remember 
that the self-field will be just compensated by the field that would obtain 
at the point O in the cavity produced by removing from the medium the 
doublet located at O. 

Now the effect at O of all the dipoles including the one at O will be 
equivalent to the field at O due to the charges developed on the surface 
of the medium due to polarizing it. Hence the averaged self-field of the 
doublet at O will be just the difference between the field at O due to the 
charges developed on the surface of the medium as a result of polarizing 
it, and the field at the same point due to all the surrounding dipoles 
in the medium. Since the difference between F and E is the polarization 
field, the latter field can be regarded as equal to the averaged self-field, 
but of opposite sign to it. 

Now the self-field of a dipole should depend only on the fine structure 
of the medium, and in particular should be independent of the external 
shape. One may, therefore, use any convenient shape like a sphere or 
an infinite parallel plate, and show that the self-field is equal to 


—(F Hf) P=—pe, es) 


where fP is the field at O due to all the surrounding doublets in an 
infinite medium. In the case of the alkali halides / will obviously be zero. 


§7. Tur Dispersion Formunar ror 4 Dense Mrpium 


We have specifically discussed till now the case when the only frequency 
affecting the dielectric constant is that of the principal lattice oscillation 
of the crystal. One good reason for confining ourselves to this case is 
that the polarization factor involved here is very definite, and is known to 
be equal to 47/3. But the main results obtained are applicable equally 
well to the case when the medium has more than one resonance frequency, 
and when some or all of these are electronic frequencies, instead of lattice 
frequencies. Let us denote these frequencies and the physical constants 

“associated with them by the subscripts Chip ke, Deheet e The polarization 
factor p,, defining the polarization field acting on any oscillator 7 due 
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to the dipole moments of all the oscillators j in the medium will not in 
general be equal to 47/3 but will depend on both i and j. The dielectric 
constant of the medium will then be given by 


e—l=4ry=—4n2 xj, - + = + ee (23) 
i 
where x; is the contribution to the susceptibility, i.e. to the polarization 


per unit volume per unit field in the medium, from the oscillators of 
type 2. 


Gs 
Amy re 1 Ea rd 
Xi Q2—a : iiX is 
in which the summation extends over all values of j namely 2, j, . . . ”. 
4an.f.e? 
G = ee 
Bi 


where n, is the number of oscillators of type 7 present per unit volume 
of the medium, /; is the oscillator strength, and y; is the effective mass 
of the oscillator. , is the frequency of the oscillator in the free state, 
to be more precise, the frequency which it will have if the influence of 
the interactions of the surrounding dipoles is eliminated. 

In the special case when all the p,;s have the same value, say p, 
eqns. (23) and (24) reduce to 


e—l B; ‘ 
— — oF 
meter bo 
where 
pC; Aor 
—— —* —— es ») 
u de p 1h . . . . . . (27) 


When further all the polarization factors have just the Lorentz value, 
namely p=47/3, it will be seen from (27) that « reduces to 2, and eqn. (26) 
to the well-known Lorentz formula. 


Herzfeld and Wolf (1925) have shown that the simple Lorentz formula 


e—l1 B; 
is mathematically convertible into the Drude formula 
C. 
=e i 
€ 2 eee <a kee Os, eaten ME) 


and vice versa. Using the same method it can be readily shown that the 
more general formula corresponding to (23) and (24) in which the polari- 
zation factors p,; may all be different, can also be reduced to the Drude 
formula (28) or to the simple Lorentz formula, and conversely the Drude 
or the simple Lorentz formula can also be transformed to a formula of 
the type (24), with any specified values of p;5. 

It will be seen immediately, whatever may be the actual polarization 
fields that may obtain in the medium, and even if there is no polarization 
field at all, that the dispersion data can be fitted, as may be desired, into 


and the Resonance Frequencies of the Alkali Halide Crystals 27 


a Drude formula, or a simple Lorentz formula, or a generalized Lorentz 
formula of the type (23) and (24) with any postulated set of polarization 
factors p;;. The characteristic frequencies that appear in these alternative 
formulae with which the same set of observational data can be fitted, will 
naturally be different ; those that appear in the Drude formula will be 
the actual resonance frequencies w, of the medium, and in the other 
formulae the frequencies will differ from w, by amounts which are deter- 
mined by the polarization fields postulated on which the formulae are 
based. Hence dispersion data by themselves can yield information 
regarding only the resonance frequencies of the medium and _ their 
oscillator strengths, and none at all regarding the polarization field. 

It may not be out of place here to emphasise these results. Though 
the mathematical equivalence of the Drude and the simple Lorentz 
formulae, which was demonstrated long ago by Herzfeld and Wolf, is 
generally conceded, the agreement between the characteristic frequencies 
appearing in these formulae, and the observed absorption frequencies is 
sometimes invoked in order to decide between the alternative formulae 
as regards their experimental validity. On considerations of convenience, 
the Drude dispersion formula, in which the characteristic frequencies 
that appear are also the observed absorption frequencies, is naturally 
preferable to any alternative formulae, but the latter formulae too’should 
be regarded as having the same experimental validity as Drude’s. 

It should be mentioned immediately that the equivalence between the 
different alternative formulae relates to the formulation of the dielectric 
constant as a function of the frequency of the incident electromagnetic 
waves, with which we are concerned here, and does not extend to the 
expressions for the dielectric constant as a function of the density of the 
medium. Since the polarization field will depend on the density the 
latter formula will be conditioned by the actual polarization field. Even 
here, the frequencies 2; of the individual oscillators, and the polarization 
field factors p;; may also be dependent on the degree of close-packing. 
We shall postpone consideration of these aspects to a later section of 
the paper. 


§8. THE Case OF THE ALKALI HALIDES 


The case of the alkali halides is exceptional, as we have mentioned 
already. The only infra-red resonance frequency that is involved in 
the dispersion formula is that of the principal lattice oscillation of the 
crystal. Making the subscript i refer to this frequency, Pi 88 we have seen, 
has just the Lorentz value 47/3. The other frequencies, denoted by 
the subscripts j, k,..., will all be electronic and the separation of the 
charges in the corresponding dipoles will be comparable with the distance 
between the neighbouring ions. Hence the distribution of the field ue 
the neighbourhood of any of the latter dipoles will differ considerably 
from that of a point-dipole. Indeed, even the order of magnitude of 
the field may become different when the approach to the dipole is close. 
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In other words, when there is considerable overlap of the electronic 
clouds of neighbouring ions, as is the case with the neighbouring ions in 
the alkali halide crystals, the polarization factors of the type pj; 
Di » ++ Will be much less than p,;=47/3, and the factors 9;;, Pex - + + 
may be still smaller. ; 

In the paper referred to, it was shown that in the alkali halides all 
the p’s except p;; have practically zero values. This circumstance, 
namely, that only one of the p’s has a finite value, and the rest have 
zero value renders the case of the alkali halides exceptional. 

The expression for the dielectric constant of these halides will, therefore, 
be of the form 
Ey—l= ema Sree Cs 

Jj 


o 2 2 
Ww; —W 


(30) 


> 9° 
w—w* 


where the electronic resonance frequencies w,;, w,, . . . will be practically 
the same as for the individual (not isolated) atoms, namely, 2;, 2,., . 
respectively, whereas the lattice frequency w, will naturally be’ very 
different from 2,. 


In particular, in the region of very long wavelengths in which the ~ 


contribution to the dispersion from the electronic frequencies, namely 
n~—1= 2 Che eecg MC ay ee ee 


is practically independent of w, and is nearly the same as for w=0, we 
obtain 


Ge 
Eg —Ny = 3 Sm: 5 niin,’ nats) ete) dae oo eee 
WwW; ——( 0's 
or alternatively, 
En—Ng" bs C,;/3 
Ey—Ny- +3 4 2 72—w?’ Se oo, 
the relation between Q, and w; being given as before (see (12) and (13)) by 
5 0; Ep — Nye +3 
Q2=07+ 3 ifs — ° . . . . (34) 


In particular one obtains from (31) for the dielectric constant in a 
static field, the well-known expression of Born 


Ge 
pty ta. er 


t 


§ 9. Discussion oF THE AVAILABLE EXPERIMENTAL DATA 


‘i The available experimental data for the refractive indices of the alkali 
‘aoe which extend over a wide spectral range, have been discussed 

y Herzfeld and Wolf (1925), by Fuchs and Wolf (1928), and more 
recently by Ramachandran ( 1947) and by Radhakrishnan (1948) 
Attention may be drawn in particular to the following results. ; 
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(1) When the data are fitted into a dispersion formula of the Drude 
type, the characteristic frequencies that appear in the formulae are 
found to agree well with the observed resonance frequencies. Since the 
fitting is done empirically, it will be more correct to say that the available 
experimental data can be fitted satisfactorily into a Drude dispersion 
formula, using the observed resonance frequencies of the crystal as the 
characteristic ones that should appear in the formula. For details 
regarding the fitting with the experimental data reference may be made 
to the original papers. The formula in general, involves three frequencies 
in the far ultra-violet, and the principal lattice frequency in the infra-red. 

(2) Since the only infra-red frequency involved in the formula is that 
of the principal lattice oscillation of the crystal the corresponding 
oscillator strength C’; will be that defined by (16), and given by (11). 
Hence it will be possible to compare the values of C, calculated from (11) 
with the numerical values appearing in the empirical formulae given by 
Ramachandran and Radhakrishnan. They are collected together in 
table 1 which includes all the crystals for which experimental dispersion 
formulae are available. 


Table 1 


Crystal LiF | NaF | NaCl | KCl | KBr 


C7 10-= from (11) 210 68 28 15 10 


C; x 10-*6 from dispersion 18 11 


(3) Now the measurements on dispersion usually extend at the infra- 
red end to about 22 1 for some of the crystals, and to about 9p only for 
most others. Since even at the former limit the contribution to n’—1 
from the infra-red term in the dispersion formula is very small, the 
experimental value of C; appearing in the infra-red terms is slightly 
uncertain, except where the value of the dielectric constant for either 
steady fields or for the usual radio frequencies has also been utilized in 
determining the constants of the dispersion formula. Actually this has 
been done only for a few crystals, by Radhakrishnan. Hence it will bs 
of interest to compare the calculated and observed values of ky—n% 
in which k, is the dielectric constant for static fields, and Muy" is the square 
of the refractive index extrapolated to w—> 0. The experimental value 
of n 2 is readily obtained from the dispersion formula by taking the 
contributions at w=O0 from all the terms except the infra-red oe The 
experimental values of ky—m9° thus calculated are given in table 2 along 
with the values calculated from (34), in which we have used for the 
frequency of the principal lattice oscillation of the cryeual the Br 
calculated by us in an earlier paper, on the basis of Born’s model, from 
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the known Coulomb and repulsion interactions between the ions_ 
(Krishnan and Roy 1951). The value as was shown in the paper agrees 
well with the observed reststrahlen frequency. 


Table 2 


ky —ny?=Ne?/zpv? 
Crystals 


Observed | Calculated 


Lif . . 7-4 
Nak . . 4-3 
NaCl : 3 
KCl : . 5 
KBr : 2: “4 


$10. DEPENDENCE ON DENSITY 


Extensive data are also available for the temperature variation of the 
refractive index of the alkali halides, and for some of them for the variation 
of the refractive index with the density also. One may therefore compare 
the observed values of (dn/dt) and —p(dn/dp)x where « is the coefficient 
of thermal expansion, in order to find whether the change of refractive 
index with temperature is due wholly to the change of density accompany- 
ing the temperature change, or whether there is in addition a pure 
temperature effect due to the thermal agitations of the ions. From the 
values of dn/dt and of —ap(dn/dp) given in table 3 it will be seen that the 
pure temperature effect, if any, is practically negligible. 


Table 3 


10°(dn/dt) 
Crystals 105 p(On/Op)_ | ——_ 
Observed | Calculated 


NaCl 12-0 0-24 —3°8 —2°9 
KCl 11-4 0-23 —3-6 —2-6 
KBr 12-0 0-35 —3-6 —4-2 
Isa 13-5 0-43 —5-0 125.8 


Now the available data for the temperature variation of the refractive 
index which extend from the extreme quartz ultra-violet to the infra-red 
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in most crystals, have been analysed by Ramachandran and by Radha- 
krishnan. They find that the temperature variation can be expressed 
empirically in the form* 


dn 2C yb, 20 ab; D 
Sh ee PP ee a es Que 
Se  ek=at* orate tara ~ 89 
where 
dw; 
a ree eee Ma ae cha ie (3) 


The first term on the right hand side of (36) is due directly to the change 
in density, since the expression for n?—1 will contain a factor proportional 
to p, and this term accounts for the bulk of the temperature variation. 
Attention may be drawn here to the extra term involving the infra-red 
frequency and proportional to 1/(w?—a?). 

We have referred to the experimental finding that in the alkali halides 
all the polarization factors are zero except p,;, which is equal to 47/3. 
Any change of density is not likely to affect p,;. If the other factors also 
remain unchanged, i.e. continue to be zero, it will be difficult to explain 
either the change in the w,’s suggested by the analysis of the experimental 
data, or the extra infra-red term proportional to 1/(w?—w?) which cor- 
responds to a change of the oscillator strength with temperature. The 
obvious inference is that though all the p’s except p;; are quite small, 
their variations with density, and therefore also with temperature, are 
considerable. Further the form of the empirical formula suggests that 
coefficients of the type (dp,,/0t), (0p;,/0t) . . . relating to the interactions 
of the electronic oscillators of the same frequency, and terms like dp,,/0t 
which define the interaction of the infra-red frequency with the other 
frequencies, are the predominant ones. These differential coefficients 
can. be calculated from the experimental data. We shall merely mention 
that dp/dt is of the order of 10~4. 


$11. FurTHER REMARKS ON THE ELECTRONIC OVERLAP OF THE 
NEIGHBOURING Ions 


In the previous sections the effect of the electronic overlap between 
the neighbouring ions has been regarded as affecting the polarization 
field, and rendering it almost nothing when the ions are at their equili- 
brium positions. Alternatively, the effect of the overlap can also be 
regarded as reducing the effective charge on an ion from e to fe where f 
is a factor which is determined by the degree of overlap, and is less than 
unity, while the polarization field is allowed to have the Lorentz value 
(Frohlich 1949). The former point of view, in which the overlap is regarded 
as reducing the polarization field to nothing, has certain advantages over 
the alternative point of view, and has therefore been adopted here. 


es ee 
* Actually these authors use wave-lengths instead of frequencies, but the 
latter are used here for convenience. 
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That in all the alkali halides the electronic polarization field is found 
to be practically zero when the ions are at their equilibrium positions may 
be theoretically significant. According to the simple Born model on which 
the foregoing discussions are based, the ions in the lattice are held in their 
equilibrium positions by the Coulomb electrostatic interactions between the 
charges carried by the ions, and the repulsion interactions between the 
neighbouring ions. When the ions have taken their equilibrium positions 
under the corresponding attractive and repulsive forces respectively, 
there will be a certain optimum overlap between the neighbouring ions, 
and this overlap is apparently also the one that corresponds to zero 
polarization field. This view receives strong support from the obser- 
vation referred to in the previous section, namely this: though the polari- 
zation field is zero for the equilibrium positions of the ions, it varies 
rapidly about this value with change of density. 
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SUMMARY 


The atomic beam magnetic resonance method has been used to find the 
nuclear spins and the ground state hyperfine structure splittings of the 
radioactive alkali atoms 74Na, “K, 8°Rb, 410s and !4Cs. The spins and 
Av values were obtained by observing 4F=0, 4dm=-+1 transitions in 
weak and intermediate fields, and the Fermi-Segre relation was used to 
calculate the magnetic moments from the results. More accurate dv 
values were also found for 24Na and 422K by locating 4F=+1, dm=+1 
transitions. The results are tabulated below. 


Isotope df H.F.S. (Mc/s) wu (NM) 
4Na 4 1139-:35+ 0-10 + 1-688 --0-005 
SAAS 2 1258-9 + 0-1 —1-137-+0-005 
86Rb 2 3960 + 20 —1-69 +0-01 
sO 5/2 13200 +110 +3-48 +0-04 
nes 4 10440 + 30 2:95 -+40-01 


§1. INTRODUCTION 


THE experiments performed by. Davies, Nagle and Zacharias (1949) using 
22Na, 4°K, 125Cs and ®7Cs, showed that the atomic beam magnetic resonance 
method is suitable for the measurement of nuclear spins and ground state 
hyperfine structure splittings of alkali atoms when only small quantities 
in samples of low concentration are available. This paper contains an 
account of experiments in which *4Na, 2k, 8*Rb, 10s and méCs were 
used in a similar apparatus. The *°Rb, #1Cs and Cs beam intensities 
were measured with a mass-spectrometer detector, and the **Na and ?K 
beam intensities were found by collecting the ionized atoms on a target 
which was then removed from the apparatus to find the activity deposited. 


§2. THE EXPERIMENTAL PROCEDURE 
2.1. The Method 


The atomic beam magnetic resonance method developed by Rabi and 
his collaborators (Zacharias 1942) makes use of the fact that the nuclear 
spin and the ground state hyperfine structure splitting of an alkali atom 
i ae ee ee aeariic unatinpaiee senor 

* Now at Glasgow University. r 
i Pe noted by Professor O. R. Frisch, F.R.S. 
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can be deduced from the way the energy difference between a pair of 
magnetic levels varies with the magnitude of the external magnetic 
field H. This is found experimentally by inducing transitions between 
the levels by applying an oscillating magnetic field perpendicular to H ; 
the energy change associated with the transition (or ‘flop’) is then 
h times the frequency at which the transition probability is a maximum. 

The experimental arrangement given in fig. 1 is used to detect the 
flopped atoms. The beam of neutral atoms, defined by the oven and 
collimating slits, passes through three magnetic fields as shown. The 
intense inhomogeneous deflecting fields A and B produce large equal 
deflections in the plane of the detector, while the field C, on which the 
oscillating field is superimposed by the RF loop, is homogeneous and 
produces no additional deflection. The three fields are in the same direc- 
tion and there is no region where H is zero or varies rapidly, so the 


Fig. 1 
2" PUMP & TRAP 0.020" WIDE TUNGSTEN 
OVEN GATE- VALVE To 2" Pump FILAMENT 
COLLIMATING SLIT RF LOOP FIXED \ 
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0 
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A MAGNET C MAGNET 8 MAGNET 


J-—_ 
2° PUMP & TRAP TO 4* PUMP REMOVABLE TARGET 


6" 
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A schematic diagram of the atomic beam magnetic resonance apparatus with 
the *4Na detector in position. - 


atoms move adiabatically from one magnet to the next. Under these 
conditions, the deflections due to the A and B magnets add up and the 
unflopped beam is split into two Stern—Gerlach components which fall on 
either side of the detector. When the oscillating field induces transitions 
between states associated with opposite effective magnetic moments in 
intense fields, the B magnet deflections oppose those due to the A magnet 
and the flopped atoms are refocused on the detector. The measured 
. beam intensity is, therefore, a maximum when the frequency of the 
oscillating field satisfies the resonance condition for a transition between 
states with opposite intense field effective magnetic moments providing 
the transition is allowed by the selection rules. 


2.2. The Energy Levels 


It is convenient to describe the energy levels of an atom in a magnetic 
field in terms of the total angular momentum quantum number F and the 
resolved component m in the direction of the field. When the total 
electronic angular momentum is 1/2, as it is for the *Sij2 ground state of 


Magnetic Moments of 24Na, 2K, 86Rb, 1310s and 14Cs 35 


an alkali atom, F takes on the two values (I +4), and the way in which 


the magnetic level W,,,, varies with magnetic field is given by the Breit— 
Rabi (1931) formula, 


hAp 4ma 12 

W m= — 2(2T-f1) +g pio +h Av (1+ eee +2) ape) 
where hAv is the hyperfine structure splitting in zero field, g;—= —p,/I, 
Hz is the nuclear magnetic moment, uy is the Bohr magneton, 
L=(97—Gr)beoH/hAv and g,; is the Lande factor for the *Sij2 electronic 
ground state. The positive sign refers to the magnetic levels associated 
with the F = (I+ 4) zero field level and the negative sign to those associated 
with F=(I[—3). The form of the level scheme for an alkali atom with 
I=4 and p; positive is shown in fig. 2; the diagram is inverted when the 
magnetic moment is negative. 


Fig. 2 


TRANSITION USED 


° ze) 20 3-0 
The energy levels of J =1/2, /=4 atom ina magnetic field. The parameter z is 
equal to (97 —91)#oH|(hAv) and the magnetic moment has been assumed to 

be positive ; the diagram is inverted when yp; is negative. 


2.3. The Measurement of Nuclear Spin 

Since the effective magnetic moment is equal to —dW/dH and the 
transitions are subject to the selection rules Ar=0, +1, dm=-+1 when 
they are induced by an oscillating field perpendicular to H, the only low 
frequency transition which satisfies both the arecHor rules and the 
refocusing condition, is the one from = (I+4), m=—(I+4) to F=(I+39), 
m=—(I+4)+1, shown in fig. 2. ; 

D2 
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The frequency associated with this transition, if H is small (7<1) and 
terms involving dv and jz are neglected, is, 
p= 2u HH /[(2I + 1)h]=2-8H]/(2I+ 1) Me/s. Fis te ge ean 
When v is less than 10 Me/s, the »z and Ay terms are small compared 
with the resonance half-width given by the uncertainty relation 
V4,2~1/t~0-1 Me/s, where t, the time spent by the atom in the RF field, 
is about 10-5 seconds. The resonance frequency is seen to depend only 
on H and J, and the frequencies associated with two isotopes in the same 
magnetic field are related by, 
v,(21,+1)=v,(21,+1). Fh ee eee 
The frequencies associated with the different spin values differ by several 
resonance half-widths when v is a few megacycles per second, so the 
method is capable of giving an unambiguous value for the nuclear spin. 


2.4. The Measurement of Av and pz 
When the magnetic field is increased until the terms involving 4v and 
g, are larger than the resonance half-width, the magnitudes of both dv 
and jz; can be deduced. If the resonance occurs at frequency v in a 
magnetic field H, Av is given by the following relation which was derived 
from eqn. (1), 

y (v—1-4g9,H) (2:800H —v) 
"~~ <2-800H 2:SigyH 
f= oT eigen ot ) 

The value of H is found by observing the resonance associated with an 
isotope of known J, «; and Ay in the same magnetic field H, and g; is 
calculated from Ay using the following form of the Fermi—Segre (1930) 
relation, 


Méjgea > 0 = aa 


Ay, 93(21,+1) 

Avg g(a 2 ee. Gale” 5. Sa 
which connects the values of J, gy and Av for two isotopes of the same 
element. The resonance frequency is found at several values of H and 
the sign of g; which gives a constant set of Jv values is assumed to be the 
correct one. 

Some accurate independent measurements of g; and Av of stable isotopes 
have shown that eqn. (3), which is based on the assumption that the 
nuclear dipole is concentrated at a point, gives g,/gy ratios which can be 
in error by as much as 0-5% (Kisinger 1952). This hyperfine structure 
anomaly does not effect the determination of the sign of g;, however, since 
the terms in eqn. (2) are small, but it does set a lower limit to the accuracy 
of a 4; value which is deduced from a measurement of Ap. 


2.5. The dF=+1 Transitions 
Several 4/=--1 transitions which satisfy both the refocusing condition 
and selection rules can also be observed in the apparatus. The pattern 
expected when H is small consists of 2 equally spaced lines with a fre- 
quency difference between alternate lines of 5:6H /(2I+-1) Me/s. The lower 
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(2/—1) lines have the same intensity and are centred on the frequency Ap 

while the other line, which comes at the high frequency end, is only half 
as intense. When J is integral, one line occurs at the frequency Av and is 
almost frequency independent ; this line was located for 24Na and 2K 
and plotted out in detail so that an accurate Ay value could be obtained. 
No corrections were made for the higher order terms which come into the 
exact expression for the frequency of this field independent line, since they 
are small compared with the error due to the wavemeter and the 0-1 Mc/s 
half-width of the resonance curve. 


§3. THe APPARATUS 


3.1. General 


The general arrangement of the deflecting system is shown in fig. 1. 
The oven and magnets are built into a 6 in. diameter brass tube and gate- 
valves are provided so that the oven and detector chambers can be let up 
to atmospheric pressure without effecting the vacuum in the magnet 
chamber. 

3.2. The Oven 

The oven used with *4Na and “K is made of monel metal and is large 
enough to contain a 2¢ charge of alkali metal. Radiation heating is 
arranged so that the 0-002 in. slit is at a higher temperature than the 
cavity, and the oven block, which contains the cavity and slit, can be 
removed without disturbing the heater coils and thermocouple. The 
block is replaced by one with a small cavity and a long channel, similar to 
that used by Davies (1949), when working with °*Rb, 1°1Cs and ™Cs. 


3.3. The *4Na and #”K Detector 


After being collimated by an 0-010 in. slit mounted on the C magnet, the 
beam falls on a 0-025 in. wide hot oxidized tungsten strip, one end of 
which is grounded. The positive ions are attracted to a 1 cm diameter 
brass target mounted on an insulated support which can be removed 
from the vacuum system through a gate-valve (fig. 1). The target, which 
is 100 v negative with respect to the ground, is connected to the input of 
a d.c. amplifier and current integrator. A second similarly situated, 
normally grounded collector is mounted in the detector chamber ; it can 
be switched to the d.c. amplifier input in place of the target which is then 
grounded. This collector is used when setting the C field and making 
other adjustments, to avoid depositing unwanted activity on the target. 
The ratio of the activity collected on the removable target during a run at a 
particular frequency to the total beam current integrated over the same 
period is used as a measure of the active isotope beam intensity. Since the 
time taken to plot out a resonance curve by this method is comparable 
with the half-lives of 24Na and #2K, allowance for decay in the source is 
made by monitoring the counter with a standard made from the material 


in the oven. 


38 E. H. Bellamy and K. F. Smith on the Nuclear Spins and 


3.4. The Mass-Spectrometer Detector 

The mass-spectrometer detector is shown diagrammatically in fig. 3. 
The 1 cm high beam enters the liquid-air cooled box and falls obliquely on 
a hot 0:020in. wide tungsten strip kept at a potential which can be 
varied from +2 500 v to +4 000 v with respect to ground. The positive 
ions which leave the strip are accelerated in the uniform field and pass 
through the two slits into the 6 in. radius, 90° magnetic field, which 
focuses the ions, according to mass, in the neighbourhood of the 0-010 in. 
final slit. The ions which pass through the slit are accelerated through a 
further 4 000 v on to the first dynode of an electron multiplier detector, 
connected to the usual amplifier and discriminator. The adjustable slit 
in front of the final slit considerably improves the resolution by defining 
the angle at which the ions are allowed to enter the final slit. An enrich- 
ment factor of 4 000 per mass number is obtained at the 133 setting of the 
spectrometer. The change over from one isotope to another is done by 
altering the voltage between the tungsten strip and ground. 


Fig. 3 
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The mass-spectrometer detector. 


A switching circuit is arranged so that a cycle of RF-on for 15 seconds 
and then RF-off for 15 seconds is repeated automatically. The RF-on 
and RF-off counts are added in separate scalers and the resonance curves 
are plotted out by finding the variation of RF-on/off ratio with frequency. 


3.5. The Deflecting System 
The A, B and C magnets have soft iron polepieces and yokes energized 
by copper tube windings insulated with degreased glass fibre yee, The 
. ae B magnets have the usual cylindrical polepieces (Rabi 1934) and 
or the A magnet, the distance between the wires of the equivalent two 
wire magnet is 0-340 in., the male polepiece radius is 0-170 in., the female 
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polepiece radius is 0-186 in., the gap is 0-081 in., the undeflected beam 
passes 0-025 in. from the male poleface and the length of the polepiece is 
4:0in. The B magnet dimensions are all twice as large. The remaining 
dimensions of the apparatus are as follows : from the oven to the A magnet 
is 8-0 in., the C magnet is 5-0 in. long with a 0-250 in, gap and the detector 
surface is 12-0 in. from the end of the B magnet. 


3.6. The Oscillating Field 

Although the C magnet length is about 12-5 cm, an oscillating field 
which extends over the central 0-5 cm is used. The resonance curve 
half-widths is then about 0-15 Mc/s, which is smaller than the separation 
between the possible low field resonances due to different spins, but 
large enough for the effects of A and B magnet current drifts to be small 
compared with the resonance half-widths, providing the C field is checked 
every few minutes. The RF used to set the C field is obtained from a 
crystal controlled 10 Me/s oscillator followed by doubler stages and a 
standard signal generator provides the variable frequency. A thermo- 
couple near the loop is used to measure the RF current which is always 
set to the optimum value found when setting the C field. The frequencies 
are measured to about 1 part in 10+ with a hetrodyne wavemeter. 


§4. RESULTS 
4.1. *4Na Results 


The oven was loaded with a one gram sample of sodium metal which 
had been in a slow neutron flux of 5x 101! n/cm?/sec for one week. The 
24Na/?®Na ratio in the source was then about 210-8 and the specific 
activity 200 mc/g. The oven end of the apparatus was shielded by lead 
to reduce the radiation hazard. A refocused **Na beam with an RF 
(on/off) ratio of between 5 and 6 was obtained. Under these conditions, 
the target activity after collecting the beam for 20 to 30 minutes was 
generally about 50 counts/minute ; this rose to 150 counts/minute when 
the RF was set to a 24Na resonance frequency. The results of the first 
successful run are shown in fig. 4, where the positions expected for 
resonances associated with spins of 3, 4 and 5 are marked, The region 
between 1-3 Mc/s and 1-4 Mc/s was studied in more detail later and the 
peak was located at 1-33 Mc/s, the position expected for a nuclear spin of 4. 

The 24Na resonances in magnetic fields associated with **Na resonances 
at 30-00 Mc/s, 60-00 Mc/s, 100-00 Mc/s and 120-00 Me/s were also found at 
the frequencies tabulated below, together with the calculated Ay values 
associated with both positive and negative g; ; the frequencies are in 
megacycles per second and the magnetic fields in gauss. 

The Ap values associated with negative g; are more constant than the 
others, so the magnetic moment is probably positive. Since the Av value 
of 1139 Mc/s came within the range of an available oscillator, an attempt 
was made to find 4F—-+1 transitions and a more accurate result for Av. 
The C field was set so that the expected spacing between the adjacent low 
field 4F=+1 peaks was 0-5 Mc/s and the frequency was reduced in steps 
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of about 0-1 Me/s, starting at 1142-0 Mc/s, until a resonance was obtained 
at 1141-40 Mc/s. Other resonances were found at intervals of 0-5 Mc/s and 
the one which was field independent was plotted out in detail (fig. 5). 
The Av value obtained from this run is 1139-35+-0-10 Me/s; the magnetic 
moment calculated using the Fermi-Segre formula with values of 
Aveg—=1771-61 Mc/s and pyg=+2-215 NM is then +1-688-+0-005 NM. 
The error given is large enough to include the effect of the hyperfine 
structure anomaly. 
Fig. 4 


TARGET 
ACTIVITY 


H = 4-28 GAUSS 


SPIN 4 


t FREQUENCY MC/S.——> T 


The *4Na low frequency resonance in a magnetic field associated with a Na 
resonance at 3-00 Mc/s. 


Table 1. The *4Na Results 


Vog H Vo4g Av(+-ve gr) Av(—ve gr) 


30-00 13-93 13-93 +0-02 1142423 1129-418 | 


60-00 78-12 29-12 40-03 11484 7 1142-23 
100-00 123-23 51-60 --0-02 11484 2 113025 2 
120-00 144-13 63-82 +0-03 11434 1139+ 2 


4.2. “K Results 
A one gram sample of potassium which had been in a slow neutron flux 
of 5X 10'* n/em?/sec for one week was used in the oven. The boo SS Ad a 
ratio in the source was then about 210-9 and the Neo: activit 
20 me/g. The procedure was the same as with 24Na: since the Sone? 
rates were down by about a factor of ten, a scintillation counter was used 


to measure the activit Bee da ' 
ees ivity deposited on the targets. The results are shown 


Magnetic Moments of *4Na, “K, ®Rb, 1310s and 134Cg 41 


The four 4F=+ 1, dm= +1 peaks were also located and the frequency 
of the field independent line was found to be 1258-9-+-0-1 Mc/s. This is in 
agreement with the Av values associated with positive gz so the magnetic 
moment of 4-K (calculated assuming Ay4.= 1258-9 Mc/s, Avgg=461-75 Mc/s 
and H3g—= 9-391 NM) is —1-137+0-005. The error given is large enough to 
include the hyperfine structure anomaly effect. 


Fig. 5 


ACTIVITY 


TARGET 


FREQUENCY MC/S ——> 


39:25 139-30 1139-35 39:40 139-45 


The field independent AF=-+1 *4Na resonance in a weak magnetic field. 
Table 2. The #“K Results 


Va5 Av(+ve gr) Av(—ve gr) 


3:00 Consistent} with /=2 
20-00 1255 +38 1233-436 
49-00 : 1262-+10 1239-+10 

120-00 1260+ 3 1254+ 3 


4.3. 8®Rb Results 


This isotope was prepared by the *Sr (da) 8*Rb reaction using 15 Mev 
deuterons on a thick target. After a 2500yamp hour bombardment, 
about 0-5 me of 8*Rb was separated chemically, using 1 mg of KCI as 
carrier and placed in the oven with 1 »g of natural RbCl. This was added 
so that the mass spectrometer and C field could be set using * Rb and “Rb 
beams; the 8*Rb/**Rb expected in the source was then about 5x 10-3. 
The rubidium beam was produced by displacing the rubidium atoms with 
sodium as described by Davies (1949). 

The low frequency resonance shown in fig. 6 was found at 2-36 Me/s 
when the ®7Rb resonance was at 3:00 Mc/s, so the spin of 8°Rb is 2. The 
frequencies of the higher field resonances located are tabulated below. 
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The Ap values associated with positive g; are consistent so the magnetic 
moment of 8*Rb (calculated assuming Avg,=6834-1 Me/s, 4gz=2°733 NM 
and Av—,=3950-+20 Mc/s) is —1-69-40-01 NM. 


Fig. 6 


BEAM 
INTENSITY 


H = 4:32 GAUSS 


22 5on 2-3 2-4 25 
The **Rb low frequency resonance in a magnetic field associated with a 8’Rb 
resonance at 2-96 Mc/s. 


Table 3. The 8*Rb Results 


Var H Vag Av(+ve gz) Av(—ve gr) 


20-00 28°37 16-13-+0-02 4330+350 | 3990+330 


40-00 56-25 32-55-0-02 3980+ 70 | 3770-4 70 
60-00 83-66 49-20 0-02 3950+ 30 | 3820+ 30 
120-00 163-21 100-51 +0-04 3970+ 20 | 3900+ 20 


4.4, 13103 Results 

A few milligram sample of !1Csg, prepared by irradiating 100g of 
barium nitrate for three weeks in a slow neutron flux of 5x 101 n/cm?/sec, 
was extracted chemically using 5 mg of KCl as carrier. This was placed 
in the long channel oven with 0-1 mg of natural CsCl and the sodium used 
to produce the beam, giving a 81Cs/183Os ratio in the source of about 
oe Ls, 

Since the #3Cs, 195s and 187Qg spins are all 7/2, a search was made for a 
significant change in RF (on/off) ratio when the oscillating field was set 
to the frequency corresponding to [=7/2 and the mass spectrometer was 
varied through the 131 setting. An RF (on/off) ratio of 4-0-0-5 was 
found, but it was almost constant from the 130 to the 132 setting and then 
began to rise steadily as the 133 position was approached, so the increase 
was contributed to Cg atoms which had been scattered into the 131 
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position. At the frequency corresponding to [=5/2, however, an RF 
(on/off) ratio, which fell away to unity on either side of the 131 position, 
was found and contributed to the 1810s atoms in the beam. A second run 

_was done to find resonances at higher frequencies and two were found 
before the experiment was abandoned because of trouble with the electron 
multiplier. 


Table 4. The 121Cg Results 


V133 A V431 Av(-+ve gr) Av(—ve gr) 


40-00 111-07 52:78+0-02 | 16750-4450 | 13050-1250 
60-00 164-02 78-75+40:02 | 15265-4145 | 13200-4100 


The Av values associated with negative g; agree so the magnetic 
moment of Cs (calculated assuming  4y4.;=13 200-+100 Me/s, 
Avy3g=9 192-6 Mc/s and py33=-+2-558 NM) is +3-48-+0-04 NM. 

4.5. 34Cs Results 

A sample of CsCl in which the 4Cs/°°Cs ratio was about 2 10-5 was 
obtained from A.E.R.E. Harwell where the 14Cs had been made by the 
1830s (ny) ®4Cs reaction. A caesium beam was produced using 10 mg 
of CsCl in the long channel oven. 

The RF (on/off) ratio was measured as a function of mass-spectrometer 
setting near the mass number 134 position, at frequencies corresponding to 
I=0, 1, 2, 3, 4 and 5, and the only significant increase occurred at the 
I=4 frequency. Four higher frequency resonances were also located at 
the frequencies shown below. 


Table 5. The 4Cs Results 


A Av(-+ ve gr) Av{—ve gz) 


56-37 02 |13750-+1500| 10350-++850 


111-07 34-56-0-02 |12030-+ 300] 10525220 
164-20 5306-40-02 |11375-+ 125] 10440+100 
314-88 | 105-62+0-02 |10975+ 30] 10440+ 30 


The Ap values associated with negative g; are consistent so the magnetic 
moment of 14Cs (calculated assuming  4?43,=10 440+30 Me/s, 
~ Arrga=9 192-6 Mc/s and jiyg3=+2-558 NM) is +2-95+0-01 NM.* 
§5. Discussion 
5.1. The Shell Model 
It has been shown (Feenberg 1949, Talmi 1950) that the known magnetic 


moments of odd—odd nuclei are generally in agreement with those 
3 eee ee seen es 


imi 1 i foation by 
* Similar results have been reported in a private communication, 
V. Jaccarino, B. Bederson and H. H. Stoke, of Massachusetts Institute of 


Technology. 
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calculated using j-j coupling and Schmidt’s values for the magnetic 
moments of the odd—proton and odd-neutron states. This agreement is 
not so good when the odd-proton and odd-neutron are in different states, 
though it can be improved by using experimental values for the moments. 
of the proton and neutron states. The magnetic moment j; of a nucleus 
with spin I is then, 
er=pPeld pS pt 1) +L +1)—J x(J yt DI/RU +1) 2] 

+pylJ a(x +l) +1LI+1)—IpJet+l2Z+1)Jy], . . - (6) 
where J» and Jy are the total angular momentum quantum numbers and 


tp and py are the effective magnetic moments of the proton and neutron 
states. 


Table 6. A Comparison between the Observed Magnetic Moments of the 
Isotopes Investigated and the Values Calculated for some of the 
Possible Nuclear Configurations 


ue Configuration 
Isotope (observed) 
Proton | Neutron 


BBS 5/2 +3:48 dsy2 —* 


>#Na 4 +1-69 | d(5/2)8 sia ds/2 
5/2 ds5)2 

d3)2 d5/2 

ds5y2 dye 


dsje fr2* 


* 


fs)2 go/2* 


P12 


d3)2* 
hyzj2 
81/2 
dsj2 
hie 


The observed and calculated magnetic moments of some of the possible 
configurations of the nuclei investigated have been tabulated below. The 
magnetic moments under the heading (8S) were calculated assuming 
Schmidt values for Hp and py, while those under the heading «(E) were 
evaluated using values of zp and yy found experimentally (Mach 1950) 
for neighbouring odd-proton and odd-neutron nuclei in correspondin, 
states. When these had not been measured, an average value taken Pail 
the plot of magnetic moments given by Klinkenberg (1952) was used 
The odd-proton nucleus ™1Cg is also included. : re 

The configurations given in the literature by Mayer (1951) and Klinken- 
berg (1952) are marked with asterisks. It is seen that in the cases of 


M agnetic Moments of 4Na, ?K, 86Rb, 210s and 134g 45 


#*K, ®°Rb and ™4Cg, the (H) values for the states are closer to the observed 
magnetic moments than the (8) values and the agreement is sufficiently 
good to be used as evidence in favour of these configurations which were 
suggested because they are most consistent with the shell model. The 
**Na agreement, however, is not good enough to decide between the 
dj. proton state predicted by Mayer (1951) and the (d5/2)%3/ State which 
is the most likely for the odd-proton nucleus 22Na. This second con- 
figuration has a.4.(E) value of 2-18 NM which is close to the observed 23Na 


magnetic moment (2-22 NM) and has been justified theoretically (Kurath 
1950, Talmi 1952). 


5.2. Decay Schemes 


The spin values found for the isotopes suggested are all consistent with 
the decay schemes given in the literature ((Nuclear Data). The **Rb 
result is of interest since the B-decay to the ground states of 8°Sr has an 
‘a’ type first forbidden shape (Zaffarano 1948, Muether 1950) which 
predicts [=2 and even parity for **Rb, assuming that 8*Sr, being an even— 
even nucleus, has spin 0 and even parity. Since the configuration 
(f5)2, S92)2, Which gives a good magnetic moment agreement, has odd-parity, 
this prediction is verified. 
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ABSTRACT 


In an attempt to increase the rate at which cloud-chamber photographs 
of V-particle decays can be obtained a number of triggering arrangements 
incorporating BF, neutron counters have been used on the Pic-du-Midi 
(2867 m). 

The yields of V-particle decays are found to be comparable with the 
yield of a typical multiple Geiger counter system. 


§1. INTRODUCTION 


One of the main problems in obtaining cloud-chamber photographs of 
V-particle decays lies in detecting the interactions in which they are 
produced as efficiently as possible. Several teams are working at 
mountain altitudes, where the flux of nuclear interacting particles is 
considerably greater than that at sea-level, but the yield of photographs 
of V-particle decays is small. Most control systems employed in this 
work are conventional arrays of Geiger counters dependent for their 
operation on the showers of penetrating particles produced in energetic 
nuclear interactions. The most serious limitation of such arrangements 
is the difficulty of detecting interactions in which only small numbers of 
mesons are produced, because of the large background of unwanted 
p-meson knock-on and air showers. It is probable that the majority 
of the penetrating shower detectors that have been used to trigger cloud 
chambers do not respond efficiently to the interactions of particles of 
energy less than about 10 Bev. Of the interacting cosmic ray particles 
with energies greater than 1 Bev less than 2% have energies above 
10 Bev. It is therefore important to consider ways in which the lower 
energy interactions can be efficiently detected. 

Methods dependent on the properties of the evaporation star accom- 
panying a nuclear interaction would appear to be the most hopeful, since 
the size of the star is not very dependent on the energy expended in meson 
production for primary energies greater than about 1 Bev (Camerini 
et al. 1951). A cloud-chamber experiment has been reported (Salvini 
1951) in which the chamber was triggered by particles whose total 


ionization loss within a crystal of sodium iodide exceeded a few Mev. 
Se a 
D * Communicated by G. D. Rochester. 
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The ionization in these cases was due to the evaporation particles from 
single nuclear interactions. This method might be rather troublesome 
to operate in a magnetic field and will probably only be useful for detecting 
interactions produced in a scintillating material. Recently workers at 
Washington University (Fowler et al. 1950, Althaus 1951 a) have reported 
the use of neutron counters, filled with enriched BF,, in coincidence with 
Geiger counters to control the operation of a cloud chamber. The inter- 
actions occurred in lead nuclei, which give a copious yield of slow neutrons. 
The yield of pictures showing recognizable mesons and V-particles was 
promising (Althaus 1951 b, Fowler 1951) and further trials seemed 
appropriate. The Blackett magnet chamber at the Pic-du-Midi (735 
g/cm?) has been used. Two major changes have been made from the 
triggering system used in the earlier neutron experiments: a large mass of 
lead on the roof of the chamber and the requirement that more than one 
evaporation neutron be detected. Asa result of the first change, evapora- 
tion neutrons are detected from a nucleonic cascade in the lead, instead of 
from just one star. The second change effectively eliminates electron— 
photon cascades, which contribute appreciably to the single neutron 
rate. The results obtained with the neutron triggering system are com- 
pared with those obtained using a conventional Geiger counter system 
on the same chamber. 


§2. A Brizr DESCRIPTION OF THE SELECTION SYSTEMS 
(i) The Arrangement Using Neutron Counters 


Four neutron counters connected in parallel were used and the electronic 
circuits, following the design of Crouch, counted events in which a number 
of neutron pulses, n, (1 >1, 2 and 3) were detected in a 400 x-sec coincidence 
gate. The gate was opened 5 .-sec after a coincidence pulse from a simple 
Geiger counter array. The relative dispositions of the counters, the 
paraffin wax moderator, the lead absorbers and the cloud chamber are 
shown in fig. 1 (a). The division of the Pb producer into two parts, with 
the C tray between them instead of being right on top of the chamber, was 
required mechanically in order to avoid altering the present thermal 
jacket of the cloud chamber. 

Various triggering arrangements were tried, such as (A,B, . Ns), 
(C,D, . N,), (C,D, . N;). Here the number gives the minimum number of 
pulses, e.g., (A, B,. N,) means one or more pulses in the A tray coincident 
with one or more pulses in the B tray, associated with three or more neutron 
pulses in the delayed coincidence gate. 

The efficiences of neutron detection were measured with a calibrated 
Ra—a—Be source. The results averaged over the volumes in question 
and corrected for the loss of neutrons captured after termination of the 
neutron coincidence gate, were 3-7% for the top layer, 2:0% for the layer 
immediately above the chamber, and 0:8% for the tungsten-copper plate 
inside the chamber. These figures show that when multiple neutron 
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coincidences are required, the efficiency of detection of nucleonic cascades 
starting in the top layer will be much greater than for the lower layers. 


(ii) The Geiger Counter Arrangement 


For comparison the Geiger counter arrangement used most recently to 
control the cloud chamber is also described. The geometry of the system 
is shown in fig. 1 (b). Two combinations in parallel have been used, the 
first requiring the discharge of at least two counters in both trays C and D 
and the second at least one in tray C and two in each of trays D and E. 


Copper 


windings 
Ae -B 
S tee 
= 
POLE Se 
| LLL LALA 5. __ LZ, 


(a) The neutron counter control system. 
(b) The Geiger counter control system, 


Areas of the counter trays: A 1000 cm?, B 240 em2, C150 em?, D 1000 cm? 
E 800 cm2?, ; 


§3. Discussion oF RESULTS 


Data relevant to the various counter arrangements described in the 
previous section are presented in thetable. During 5 
LU eucenanG p ve vble. During the tests of the neu‘ron 
gering angements, the resetting time of the chamber was two 


minutes. The data obtained with the multiple Geiger systems have ben 
corrected to this resetting time. 
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The first neutron counting arrangement (A,B, .N;), is insensitive to 
neutron-induced events, and responds mainly to interactions produced by 
protons of several Bev energy. The geometrical efficiency of the chamber 
for the detection of the products of such interactions when they occur in 
the top producing layer is small. In order to increase the yield of events 
giving particles in the chamber, an arrangement with the Geiger telescope 
below the main producing layer was used. One energetic secondary 
particle was then required below the cloud chamber. This change 


Comparison of Various Control Systems 


(1) (2) (3) (4) (9) (6) (7) (8) 
Counter | No. of | Rate | Total | No. p.p. |(a) No. of |No.hrs.| No. p.p. 
arrange- | photo- | per | time | per hr. | V°-part- | per V-| per V- 


ment. graphs.| hr.* | (hrs.) icles. part- | part- 
(b) No. of icle. icle. 
V +-part- 
icles. 
A, B,.. Ns 587 ll 53-5 3640-2 |(a) 2 ~27 ~97 
(b) 0 
C,D,.N,; | 1946 9-1 | 214-0 6-0+0-1 |(a) 11 19+6 | 114-+38 
(6) 9 
CDi oN, |, 1328 7 13-6 97-6 | 10-0-+0°3 |(a) 6 16-7 | 160-+70 
(6) 0 
C,D,.N, | 804 6-3 | 127-5 5-0+0-2 |(a) 6 18--7 | 90-435 
(b) 1 
C.D, 
and 


U,DoKs 3923 8-0 | 490-4 6-4+0-1 |(a) 30 13+3 | 83-20 


* Resetting time of chamber was 2 min. 


raised the average energy of the events detected, and also made the 
arrangement sensitive to interactions initiated by neutrons as well as 
protons. The second, third and fourth rows of the table give the results. 
Several other neutron triggering arrangements were tried briefly, but the 
results were not considered to be of sufficient significance to tabulate. The 
last row combines the results of two runs with a penetrating shower selector 
of the conventional type, between which the neutron trials were made. 

No definite conclusions can be drawn regarding the absolute yields of 
V-particle decays, given in column (7) of the table, since the numbers 
recorded are small. The first arrangement may have the lowest efficiency 
while the other neutron arrangements give yields comparable to those 
obtained with the multiple Geiger systems. 

The yields of shower particles (p.p.) penetrating the plate across the 
chamber are shown in column (5) and the numbers of penetrating particles 
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for each V-decay in column (8). Although a significantly larger number 
of shower particles per hour are detected with the arrangement (C,D, . N,) 
there is probably no significant increase in the rate of V-particle decays in 
the chamber. Moreover there are no obvious differences in the energies 
of the V-particles recorded by the neutron counter and Geiger counter 
arrangements. It thus appears that there is no great advantage in using 
a neutron counter triggering system at mountain altitudes unless the 
chamber resetting time can be considerably reduced. However, even for 
large chambers the neutron detection system can probably increase the 
rate of recording V-particle decays at sea level. 
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SUMMARY 


Proofs are given of theorems stated in the Bishop—Hill theory of 
crystal distortion. Mathematical relations are established for the 
deformation of crystals by simultaneous shears, and the possibility of 
enforcing an arbitrary strain by glide is examined. The uniform stress 
states not exceeding the critical shear stresses in the glide-directions are 
discussed, and a geometrical picture of stress and strain states is developed. 
It is proved that it is always physically possible to enforce a given strain 
by glide provided it is geometrically possible. The uniform elongation 
of a face-centred cubic crystal is examined. 


§ 1. INTRODUCTION 


_ EXPeRIMENTAL investigation of the plastic deformation of metallic and 
ionic, crystals has led to the conclusion that two basic processes are 
involved, namely, gliding and twinning. In glide, planes of atoms slide 
relative to each other through integra] numbers of lattice spaces. This 
preserves the lattice structure, but constitutes a discontinuous shear 
at the atomic level. Mechanical twinning can be regarded as a continuous 
shear in portions of the crystal, which results in the atoms taking up a 
symmetrical position with regard to the atoms in the original lattice. 
The magnitude of the shear is such that the lattice structure is again 
preserved. Both processes are commonly observed to occur in plastically 
deformed crystals. Strictly speaking, these processes are not directly 
applicable to real crystals, which are only perfectly ordered over small 
regions. However, the departure from the idealized model is sufficiently 
small to allow of real crystals being treated as perfect. The spacing 
between the individual slip-planes, and the amount of slip occurring 
on each, is such that a uniform macroscopic shear can conveniently be 
regarded as a uniform continuous shear for the purposes of a mathematical 
theory of crystal distortion. Assuming that glide alone is a sufficient 
mechanism for deformation, it has been shown (Bishop and Hill 1951 a) 
that two extremum principles hold for single crystals. The extension 
EE 
* Communicated by the Author, 
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of these principles to a polycrystalline aggregate leads to the mathe- 
matically important result that the plastic potential and the yield 
function coincide for an aggregate (as they do for single crystals). This 
in turn enables the yield to be computed. In this paper, the circumstances 
under which glide alone is, in fact, a sufficient mechanism, are investi- 


gated. 
§ 2. GLIDE-STRAIN RELATIONS 


It is known that the atomic planes on which relative sliding occurs in 
crystals are closely related to the lattice structure. In metals containing 
close-packed planes, glide occurs on these planes, and the close-packed 
directions are glide-directions. The glide-planes and directions can be 
determined from suitable surface markings on test specimens (see, for 
example, Taylor and Elam 1925). 

If glide is regarded as a continuous shear parallel to a glide plane, the 
displacement of an arbitrary point due to a given shear can be found by 
simple geometry. Knowing the displacement components, the strains 
are immediately obtainable by differentiation. 

Considering the general problem, let there be 7 glide planes with 
normals n;“, «=1, 2,...7r; 7 takes the values 1, 2, and 3 and gives the 
direction-cosines of the normal in an orthogonal Cartesian space. Let 
the glide plane « contain s glide-directions 17°, B=1, 2,...s. (Generally 
_ 8 will be different for each plane.) The displacement components wu, at 
a point P, position x;, due to a uniform shear of magnitude dy* parallel 
to the « plane in the f direction, are 


Uj =X dy" n,Mal 2, 
| 


where dy*’ is the engineering shear equal to twice the mathematical 
shear, and the origin is chosen to be a point of zero displacement. The 
summation convention of summing over repeated indices does not apply. 
The total displacement components due to shears on all planes are 
therefore 
w= 2 dye nepal.” eet ae oe 
«Bj 


Since the strain components and displacements are related by 


__lfdu; , du; 
deo 5 (Fat Be) 


the strain components corresponding to (1) are 


de y= 2 dy (Lense) He see. LOUD) 
on 


Equation (2) represents six relations since de,;=de;; In shearing, no 
volume change occurs, and it follows. that dey,+dey.+de,,—0. An 
arbitrary strain in an incompressible material is therefore specified by 
five independent strain components. Equation (2) can be considered 


of the Plastic Deformation of Crystals by Glide 53 


as a set of linear algebraic equations expressing the strains de;, in terms 
of the shears dy”. 


If the column vectors [de], [dy] have components 


{de,, OL RUE 2déos, 2d€s1, 2de1,} and {dy}, yee Oy andy, eh, } 


? 


and [H] is the matrix of the coefficients of dy*®, eqn. (2) becoms 
[de]=[E][dy]. SS Mga ee aan ACY 


Omitting one of the equations in de;,, (de,, say), if a (55) minor of [EZ] 
has a non-zero determinant, eqns. (2) contain a linear independent set 
of five. This implies that an arbitrary strain can be imposed geo- 
metrically in terms of five independent shears, i.e. by glide alone. 
(A distinction between geometrically and physically possible sets of shears 
will be made later.) Since only five independent strain components 
exist, only five shears can be linearly independent. 

Normally it is necessary to examine [H] to find whether a non-zero 
minor exists, but some general remarks can be made. Since only two 
independent shears can occur on one plane, at least three non-parallel 
glide planes must be available. (The existence of three such planes, each 
with two independent shears does not, however, imply that a general 
strain can be imposed by glide alone.) 


§ 3. SrrEss IN A SINGLE CRYSTAL 


It is known that glide can occur when the resolved shear stress on a 
glide plane in a glide direction reaches a specific value, the critical shear 
stress. This value may differ for different glide directions and for 
opposite senses of shear in the same glide direction. In considering the 
variation of the critical value with strain, two types of glide systems 
must be differentiated, active systems on which slip is actually occurring, 
and latent systems which are not slipping. The experimental evidence 
as to how the critical values on the two types of system depend on 
strain is slight. Taylor and Elam in their work on aluminium (1925) 
have concluded that both active and latent systems harden by approxi- 
mately the same amount. More recently, Rohm and Kochendorfer (1950), 
on the basis of their experiments on aluminium crystals, have proposed 
an empirical law relating the hardening on active and latent systems 
through the lattice geometry. Any applied stress system must be such 
that the resolved shear stress in each of the glide directions does not 
exceed the current critical value. 

Let 7 be the resolved shear stress on the « plane in the positive B 
direction. If o,, is the current uniform stress state, . 


Chuan a7 of 
ie = 20 Nj l; 5 
WY 


=Sho y(n {lPt+njl”). MS raat. cee) 
7) 
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This set of equations expresses the resolved shear stresses in terms of 
the stress components. If 7,%°, 7, are the critical shear stresses In 


the positive and negative senses of shear, 
wt 71% <7 <7, 
i.e. — 7% <Zho y(n 7? +771) et ae . . . . ° (5) 
ij 


for all « and f. 

If column vectors [7] and [o] have components {711, 71, ... 771, 7, .. + } 
and {0,1, G29, 633) 23, %51» Fig}, and [S] is the matrix of coefficients of 
. the o;; in (4), 

[r7]—[ Slo]. 2 oe eo, os eee ee 


(It will be realized that the coefficient of o,,,, say, will contain the terms 
in both o,, and a,,,.) Comparison of (2) and (4) shows that the matrix 
[EZ] is the transpose of [S]. This can be predicted from a different 
viewpoint. The virtual work theorem states that the work done in a 
strain by a stress system, is equal to that done by an equivalent set of 
shear stresses in an equivalent set of shears, i.e. 

20 5 deg=—at™ dy, 

aj oB 
The stresses and strains do not necessarily correspond. In matrix 
notation this becomes 

[o]'[de]=[7] [dy], 
where the prime denotes transpose. 
Now [de]=[E][dy] and [r]=[S][o]. Hence [o]'[E][dy]=[e]’[S]'[dy]. 


This is true for all [o] and [dy] and hence 


[Z]=[S]’. 


§ 4. GEOMETRICAL REPRESENTATION 


So far, we have seen that for a general strain to be geometrically 
possible, five independent shears must exist, and that mathematically 
this requires a (55) minor of a certain matrix to have a non-zero 
determinant ; further that a set of inequalities places conditions on 
possible combined stress states which may be imposed on the crystal. 
In this section it is proposed to develop a geometrical picture illustrating 
the relations between stress, strain, shear, etc. 

Consider an orthogonal space of six dimensions with (de,,, . . . 2de,.), 
(1;, - . + O19) plotted along the co-ordinate axes. Stress and strain 
states appear as six-vectors in this space. Since the work done by a 
stress o,, in a strain de,, is 


dW =o, dey, +099 dego+o55 déx3-+ 2095 déy3-+ 205; dés1+ 2049 dezp, 


the work done is the scalar product of the stress and strain vectors. 
c F veces 

Since each glide direction can be operated in two senses, twice as many 
vectors corresponding to unit shears will exist as glide directions. Any 
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five independent vectors corresponding to incompressible strains (taken 
in both senses) will map the entire subspace defined by de,,+deg.+de3,—0. 
Hence all sets of five independent shears can constitute an arbitrary 
strain from purely geometrical considerations. Any five independent 
shears operating in one sense only will map a region of the subspace. 

Now let us consider stresses. From eqns. (2) and (4) it can be seen that 
the resolved shear stress in a glide direction is the scalar product of the 
stress vector with a vector corresponding to a unit shear in the glide 
direction. Each of the stress states which achieves the critical value 
in a glide direction will therefore lie in a hyperplane normal to the unit 
shear vector and given by either 


7 = Zo ign cL or —7,'%=Zo gel. 
aj 1) 

The set of these planes constitutes a bounding surface called the yield 
surface, within or upon which all stress states not exceeding the critical 
shear value in any glide direction must lie. If such planes meet in a 
point, the corresponding stress vector will achieve the critical shear 
stress in all n glide directions, and will make the simultaneous operation 
of all m shears possible. Since a pure hydrostatic pressure gives no 
shearing stress on any plane, the yield surface will be cylindrical with 
generators in the direction (1, 1, 1, 0,0, 0). A stress vector can always 
be found to achieve the critical value in five independent shears, but 
in general it will lie outside the yield surface. We have seen that under 
. certain conditions sufficient shears will be available to impose an 
arbitrary strain from purely geometrical considerations. It remains to 
show that a stress system satisfying the yield conditions can always be 
found to operate such a set of shears. 

Geometrically the problem appears as follows. If there are N glide 
directions there will be 2N vectors corresponding to unit shears (regarding 
positive and negative shearing along a glide direction as separate). 
Suppose several sets of five independent vectors exist. All positive 
combinations of the vectors in a set will map a region of the space, 
within which all strain states imposable by simultaneous operation of 
the shears must lie. To each set will correspond a stress vector whose 
scalar product with each vector in the set gives the (positive) critical 
value. In general, an arbitrary strain will lie in a region generated by 
several sets of independent shears. It is required to show that the stress 
vector corresponding to at least one such set lies upon the yield surface. 


§5. ExIsteNCcE THEOREM 


Suppose n sets of five simultaneous shears exist which will impose 
a given strain geometrically. Consider the corresponding set of n stress 
vectors which operate these sets of shears. The minimum shear principle 
for a single crystal (Bishop and Hill 1951 a) can be stated in the form that 
if one of the vectors terminates on the yield surface, it will do less work 
in the given strain than any of the set terminating outside the surface. 
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Consider the set of shears for which the work is a minimum, and the 
corresponding stress vector. Either this lies upon the yield surface or 
else the critical shear stress is exceeded in some other glide direction. 
Let the set of (unit) shears be the vectors dy, («=1, 2,...5). Since 
we are treating positive and negative directions of shear as separate 
vectors, the given strain de can be expressed by 


de= La, dy,, 


where all the a, are non-negative. 
If +, are the critical shear stresses and o is the stress state, 
Go. dy,= ls 
Further if dy, is the direction in which the critical value is exceeded, 
o. dy, >7.. 


It is proposed to show that a positive combination of this shear and 
four of the original five will suffice to constitute the given strain, i.e. 


de=2a, dy 2=242.4Yp) G4 = sel ee 
C4 B 


where B=(1, 2,...6) with one of the indices (1, 2,...5) excluded, and 
all ag>0. 
Since the shears dy, are independent, 


dy,=2b, dy,, a=1,2,,..93 
a 


“the b,, may be positive, negative or zero. 
Hence dy ,=1/b,,{dy,—2(b, dy,,);.}; where k takes the values (1, 2, .. . 5) 
a 
and the subscript k outside the bracket indicates that k is omitted in 
the summation. 
a 
matte dy,.=2(a,, dy,,) Pa ie {dy,—2(b,. dy.) i} 


“« a 


_% My Oy 
= Fytteh ze Beg 


Since the stress vector has a positive scalar product with all six unit 
shears, and dy, is a positive linear combination of the other five, not 
all b, are negative. Consider the set b,/a,, «=1, 2,...5. Since all the 
a, are positive, at least one quotient will be positive and there will be a 
maximum value b,/a; (not necessarily unique). 

Taking k=j, it is seen that b,(a,/b,—a,/b,) is always non-negative, 
and hence the given strain can be imposed by a positive combination 
of dy, and four of the original five shears. Let 6, be the stress required 
to achieve the critical shear stresses in this new set of shears. From (7) 
and the virtual work theorem, it follows that the work done by o is the 
same in both the new and the original sets of shears. Since 6 exceeded 
the critical shear stress in dy, and o,, just attains it, 


o, -de<oca.de, 
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i.e. less work is required in the new set of shears to impose the given 
strain. But by hypothesis the original set corresponded to minimum 
work, and hence the assumption that o is not a yield state must be 
wrong. The following theorem may be therefore enunciated. If a 
single set of five independent shears exist amongst the available glide 
systems, any arbitrary strain (with zero hydrostatic part) can be imposed 
both geometrically and physically. 


§6. YIELD SuRFACE OF A CRYSTAL 


The maximum work principle for a single crystal states that the actual 
stress, corresponding to a given strain does not do less work than any 
other stress state satisfying the yield conditions. Thus to find the stress 
required to operate a given strain, it is only necessary to find the stress 
vector lying on the yield surface maximizing the work. The yield surface 
is defined by an infinity of vectors, but it is not necessary to consider 
the infinity in maximizing the work. As we have seen, the stress surface 
is bounded by plane faces. If the strain vector corresponds to a single 
shear, any stress vector lying in the appropriate plane face of the yield 
surface will maximize the work, i.e. will operate the shear. A general 
strain requires five shears, which will be operated by stress vectors lying 
at the intersection of five plane faces. Since such stresses can also operate 
any selection of the relevant shears, it is clearly sufficient to consider 
only the ‘vertices’ of the stress surface in maximizing the work (as 
opposed to satisfying the stress equilibrium equations). The term 
‘ vertex ’ is used by analogy with polyhedral surfaces in three dimensional 
space, but care must be exercised in identifying the geometrical properties 
of the vertices in the two spaces too closely. In the discussion below, 
the stress states defining the vertices of the yield surface of a face-centred 
cubic metal are derived. 

In a face-centred cubic lattice, four glide-planes with three shear 
directions in each exist. The resolved shear stresses in terms of the 
stress components referred to the cubic axes are (multiplied by a factor 


1/6) 


A-—G+H (a,) B+F—H (a,) C—F-+G (as) 
A+G-+-H (6;) Bi (be) C-+F—G (bs) 
A+G—H (c,) B+F+EH (¢) C—F_—G (¢s) 
ASG Hidy) B—F-+H (d,) O+F+G (ds) 


where A=01,—033, B=o33—011, C=041—022, P=023, G=o531, H=03,, 
The letters and figures in brackets refer to the plane and glide direction 
respectively. (See Bishop and Hill 1951 b.) F: 
In accordance with observation, we shall assume that the critical 
shear stress is the same on all glide planes, and for convenience it will 
be taken as 1/,/6 in both senses of shear. The yield surface is therefore 
defined by 
Wee Hel, BLP TH=-+1,” CEF2GSa1, 23) 


where all combinations of signs are possible, giving 24 equations in all. 
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A stress state not exceeding the critical shear stress in any glide direction 
will satisfy the inequalities 


|ALG+H|<1, |BEFLH|<1, |CHF+G|<1. 


Consider three rectangular Cartesian spaces with axes (A, G, #), 
(B, F, H) and (C, F,G). The stress components must be such that the 
vectors with components (A, G, H), (B, fF, H) and (C, F, G) lie within 
or upon the octahedra given by (8). One such octahedron is shown 
in fig. 1. The faces are labelled to show the shears corresponding to 
stress vectors terminating in them. Vectors terminating at P,Q, R 
operate 1, 2, and 4 planes of type 1 respectively. Suppose the (A, G, H) 
vector terminates at P. This fixes the values of G and H, and it follows 
that the vectors in the other octahedra lie in sections by the places 


Fig. 1 


H 


G=const., H= const. These are shown in figs. 2(a) and 2(b) for 
A=1—g—h, G=g, H=h. Since A +B-+C=0, these two figures may be 
superposed as in fig. 3. A stress state not exceeding the yield conditions 
will be such that the vector in fig. 3 lies within the common shaded region. 
The sides are labelled to show the corresponding shears operated for 
g,h>0. When g, h=0, the shears in parentheses also operate. 

If (l—g—h), g, h lie within the open interval (0, 1), only one plane of 
type 1 operates. From fig. 3 it can be seen that when g=} or h=+t 
four more shears can be operated. Since, however, these are medal 
in the sets operated by g=} and h=} (ie. A=0), the stress states do not 
represent true vertices. It is therefore only necessary to consider stress 
states with at least one of A, G or A zero. 


(i) Let A=0. A maximum number of four planes of type 1 operate. 
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The corresponding sections of the other two octahedra with g=0 are 
shown in fig. 4 (a). All possible cases are covered by g=0, g=% and 
g=1. The corresponding stress states (including the synonymous states 
in B=0 and C=0) are Nos. 4-6 and 13-28 inclusive in the table. 


Fig. 2 
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Fig. 3 


(ii) Let G=0. A maximum number of four planes of type 1 operate. 
Figure 4(b) shows the corresponding sections of the other octahedra 
for h>0. All possible cases are covered by h=0, h=} and h=1, giving 
with synonymous states for H=0, Nos. 1-3 and 7-12. 
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This set of stress states (together with the negative values) defines 
the vertices of the yield surface. They were first determined by Hill 
(1949)* by means of a different method. Each of these states can operate 
either six or eight simultaneous shears, and in general several combinations 
of these will impose a strain both geometrically and physically as will 
be shown later. The yield surface for any lattice can, in principle, be 
similarly determined. If differential hardening occurs with strain, the 
yield surface will change in both size and shape. 


Fig. 4 
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§7. Untrorm ELONGATION IN A F.C.C. Crystan 


As an example of the use of the maximum work principle in conjunction 
with the tabulated stress states, it is proposed to examine the elongation 
of a single crystal with the transverse strains constrained to be equal. 
All orientations of the crystal cubic axes relative to the specimen axes are 
treated. This problem has been considered by Taylor (1938 a and b), 
who solved the sets of simultaneous equations expressing the strain in 
terms of five independent shears, and selected as the actual shears the 
set having a minimal sum of the absolute shear magnitudes. 

If (x,y,z), (#, y’, 2’) are the cubic and specimen axes respectively, 
and x’ is the direction elongation, 


de 0 0 
de'j=| 0 —dde 0 
0 0 —tde 


* Private communication. 
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y’ and 2’ are disposable orthogonal directions normal to a’; y’ is taken to 
lie in the plane z=0. Taking 4 as the meridional angle of the plane 
containing the x’ and z axes from the plane y=0, and 7/2—¢ as the angle 
between the x’ and z axes, symmetry considerations show that all 
orientations are covered by 


0<0<7n/4, 0<¢<tan“! (sin 6). 
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If a,, is the direction-cosine of the i specimen axis relative to the j 
cube axis, | 
cos ¢ cos 6 cos ¢ sin 0 sin d 
Cee ain cos 8 0 
—sin cos 0 —sin ¢ sin @ cos 
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The strain tensor relative to the cubic axes is given by 
de p20 pies de’ 1, 
whence 
3 cos? ¢ cos? 0—1 3cos*¢sin#@cos@ 3sin ¢ cos ¢ cos 6 
de j= a 3 cos? ¢ cos@sin# 3cos*¢sin?6—1 3sindcosdsind |. 
3sin¢dcosdcosé 3sindcos¢sin#d 3sin? ¢—1 


The work done by a stress o;; in a strain de ;; is 


ij 


Since dé» +déy,+de,,=9, 
dW=—B dezztA deyy+2F de,,+2G de,,+2H dez,y. 


From the maximum work principle, the yield stress state required to 
operate a given strain maximises the work. We have seen that the tabu- 
lated stresses form a complete set for this purpose. Hence the work done 
in a given orientation (specified by @ and ¢) is the maximal value of the 
works done by the fifty-six stress states. It follows from the existence 
theorem that this state (not necessarily unique) operates sufficient shears 
to enable the strain to be enforced. 

In fig. 5 the relevant stress states are shown: the bar above the 
numbers indicates that the stress state is the negative of that in the 
table. The boundaries between states are found by equating the works 
required by the states on either side. The work required to impose the 
strain is therefore known for all 6 and ¢. In region 6 for example, the 
work required is 


dW =3 cos? ¢ sin @ cos 0 de. 
In this region, the strains and shears are related by 
/6 deyj= —A,—b.—C.—dp, 
/6 deg.=a,+b,+¢,+d,, 
24/6 des, =a,—b,+cC,—dy, 
24/6 des;=—a,+b,+¢,—d,, 
2/6 dey,=a,—a, +b, —b,—¢,+¢,—d,+dy. 
Treating a,, a, and b, as parameters, 
2¢,= 4/6 deg, +24/6 des; —2b,, 
2d,=+/6 deg.—2/6 deg,—2a,, 
2b2=— 1/6 de,,— 6 degy— 24/6 deyy+20,+ 2b,— 2a, 
2¢,= — 1/6 de,,+ 2/6 den,—2a,, 
2dy= 1/6 degs—24/6 deg, +24/6 dey.—2a,—2b,+ 2a. 
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The stress system operates the shears in the senses +@,, +b,, —c,, 


dj, —G, —b,, +¢,, 4 dy and hence a,, a, and b, must be chosen to 
satisfy these conditions. 


de iG 


Oo Be BH 
SO BR Blo 


and C= der/6/8—b,, 

d,=der/6/8—a,, 

b,=—der/6+a,+6,—de, 

Co= —der/6/8—Qp, 

dg—de71/6/8—a,—b, +4. 
Clearly a, >de1/6/8, b; >der/6/8, ay <dev/6/8, with nd tat aa 
Provided a,, a, and b, satisfy these inequalities, this sys em 0 a 
shears will impose the strain and is operable by a yield stress state. 
therefore that this particular strain may be imposed by only 


te set of shears is available. It is 
y of sets of shears 


It appears ' particule 
one stress state, but that a triply aaa pent, 

nerally true that a single stress state and a m a ) 
ae te an arbitrary strain (although certain strains may be imposed 


by several stress states). 
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ABSTRACT 


A study is made. of regular assemblies of species A and B, whose 
configurational energy is given by the expression WN 43(1)+w.N 4p(2), 
Nap(l) and N,,(2) being the numbers of first and second neighbour 
A-B pairs respectively. Positive and negative values of w, and w, are 
treated and results apply to phase-separation, order—disorder and Ising 
ferromagnets and antiferromagnets at zero field. Generalized quasi- 
chemical approximations giving the correct ratio of total first and second 
neighbour pair numbers are used and the consistency is checked by the 
critical points in the limiting cases w,=0 and w,=0 and by the value of 
w,/w, at which the critical temperature becomes zero. The simple cubic 
lattice is treated in detail and a critical temperature curve derived for the 
plane square lattice. Fluctuation effects are diminished for reinforcing 
first and second interactions and increased for opposing ones. The 
quadruplet approximation which has been used by certain authors is 
critically examined and it is pointed out that it would be consistent for the 
disordered state on the face-centred cubic lattice. 


§1. INTRODUCTION 


THIS paper will be concerned with the statistics of a regular binary 
assembly in which the interactions between first and second neighbours 
are appreciable. A method of constructing a generalized quasi-chemical 
approximation into which second neighbour interaction energies are 
incorporated will be described and certain criteria of consistency by which 
such an approximation can be tested will be developed. The assembly 
on a simple cubic lattice will be dealt with in detail. 

In real condensed phases further neighbour interactions are not 
negligible, though they have often been assumed to be so in statistical 
studies. Even with an inverse sixth power law the ratio of the energy of a 
second neighbour pair to that of a first neighbour pair is 1 : 8 for a simple 
cubic lattice and 27 : 64 for a body-centred cubic lattice. It is considered 
that agreement between the computed and observed effects of superlattice 
formation could be considerably improved by taking account of the second 
ee eee 

* Communicated by the Author. 
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as well as of the first interaction (Chang 1937, Fournet 1951). In materials 
where there are forces of the exchange type the further neighbour energies 
may be comparable with, or even larger than, the first (Seitz 1940, 
Bozorth 1947, 1952). A statistical method which gives proper weight 
to the effect of second neighbour interactions is therefore of considerable 
interest. 

Consider a lattice of N sites, each occupied by either a system of species 
A or a system of species B. Let the number of A be NV, and the number 
of B be Nz so that N,+N,=N. Let each system have z, first neigh- 
bours and z, second neighbours so that there is a total of 3z,N first 
neighbour pairs and 4z,N second neighbour pairs. Let the number of 
‘unlike’ first neighbour pairs be N,,(1) and the number of ‘ unlike’ 
second neighbour pairs be NV,;(2). The configurational energy, #, is 
linearly dependent on the various pair numbers and by a suitable choice of 
zero can be put into the form, 


B= wv Nanh)twe ales. soe. ee ee 


where 2(z,w,-+zW,) is the energy increment on exchanging an A and a B 
between pure A and B phases. 

In the generalized quasi-chemical method (Guggenheim 1944, Yang 
1945, Guggenheim and McGlashan 1951 a) an approximate configurational 
partition function is constructed by using a basic group of n sites, in which 
there are m, first neighbour pairs. To get the right total number of first 
neighbour pairs an assembly of 4Nz,/m, of these groups is used and, 
assuming their independence, a combinatorial factor is formed in terms of 
the probabilities, y,, of the various configurations, 7, on the sites of group. 
This is then multiplied by a normalizing factor g)(NV,, Vy), given by 


In 9o(N 4, Np)=($nz,/m,—1)(N, nN,+NpglnNg—NInd). . (2) 


This is determined from the partition function when the interactions are 
zero and will thus be unaffected by the introduction of second neighbour 
energies. 

Now suppose that the group contains m, second neighbour pairs ; the 
number of such pairs in $Nz,/m, groups will be }Nz,m,/m, and for this 
expression to be equal to $z,N it is necessary that 


M>/M,=2,/23. es te ee 


Consistent groups for which the relation (3) is satisfied will be used in this 
paper. Several attempts have been made (Chang 1937, Hill 1950, 
Guggenheim and McGlashan 1951b) to extend the quasi-chemical 
approximation by using groups in which m,/m,4z,/z, and these will be 
discussed below. As ‘a posteriori’ criteria of the consistency of the 
approximation the critical temperature values in the limiting cases w,=0 
and w,=0 and the value of w,/w, at which the critical temperature tends 
to zero will be used. 
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If w,=0, w,>0 and N, and Ny are fixed the model represents a 
disordered binary mixture. If the lattice divides into equivalent sub- 
lattices a and 6, such that all nearest neighbours of a given site are in 
the other sub-lattice, then, it is well known that this is equivalent to 
another binary mixture, with N,—N,, and a negative interaction —Wy. 
Let the ratio N ,/(N,+N,) be denoted by x and the degree of order with 
respect to sublattices a and b by s, $(1—s) being the proportion of systems 
in ‘wrong’ positions. Then the compositions of the conjugate phases of 
the disordered mixture are related to the equilibrium degree of order in the 
other model by 


Pea has). eo l=e—t(I—s) (4) 


This can be shown by replacing each A ona b site by B and each B by A, 
leaving the configuration on sub-lattice @ unchanged. If the lattice is 
such that all second neighbours of a given site are in the same sub-lattice 
as this site then this interchange will not affect the total second neighbour 
energy, W.N,4,(2). A model with interactions w,, w, is thus equivalent 
to one with interactions —w,, w. and eqn. (4) still applies (Bell 1952). 
Define a parameter p by p=w,/w, where w,>0 and p=—w,/w, where 
w,<0. The first and second interactions will now reinforce each other in 
both models when p>0 and oppose each other when p<0. The mathe- 
matical theory can be developed in terms of either model; the positive 
first interaction (disordered) case will be chosen here. 

When JN ,, Ny» are not fixed the disordered model becomes a zero-field 
Ising ‘ferromagnet’, which is quantum-mechanically incorrect but 
permits statistical calculation over the whole temperature range. (For a 
full discussion, see Van Vleck 1945.) The relative magnetization is 
connected with the conjugate phase compositions by relations identical 
with (4). The ordered model is related to an ‘ antiferromagnetic ’ model 
of similar type. (For the ‘ first neighbour only ’ case with applied field, 
see Brooks and Domb 1951, Ziman 1951.) 


§2. THE StimpLE Cusic Lattice 


In this lattice z,=6, z2,=12, so that the ratio mg/m, for the basic group 
must be equal to 2. The group employed consists of a central site with its 
six nearest neighbours so that n=7, m,=6, m ,=12. It will be repre- 
sented schematically by 


where ‘0’ is the central site, the ‘ X’ are its four nearest neighbours in 
one plane and ‘+’ are the other two nearest neighbours. There are 


F2 
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twenty types of configuration on these sites and the first ten of these, with 
their ‘ weights ’ w,; are given below : 


(1) A (2) B (3) A (4) B (5) A 
A B A B A 
AAA BoA 5 BAA A AB Be At 
A B A B A 
A B A B A 
o.= Wo=1 wWs—6 O,= on= 
(6) B (7) B (8) A (9) A (10) B 
B B A B B 
AAA AAA BA B Bab BAA 
A B B A A 
A B B A A 
w@g=12 wW7=3 wg=12 Wg=12 W19=8 


Configuration i+10 is obtained from i by replacing each A by B and each 
B by A. 

In order that the probability of occurrence of an A on both the central 
site and the outer sites should be zw, and that the sum of the configuration 
probabilities, Y,« /,, should be equal to unity, the following three relations 
exist. 


a= ys, Wo Bysg+ Gyhy + 8h + 12g + 37+ L2pgt l2yg+8%y9, . - (5) 
Ya yy + yn + 843+ Spat Beq5 + 12g + 817+ 1213+ 1219+ 829, (6) 


Pops t dpa tpst det 2h, + Spe t 69+ 419 
=Yyot Past Spat past spre t 212+ Spas t S19 + 429. 6 Wed tae 
The pair numbers NV, ,(1) and V,,(2) are given in terms of 4; by 
N gpl) =3N (po+ 2+ 3+ 13+ dehy Sify 4 ths +s t+ thet 4g 
+ eq + 2th, + 8ipfgt Siphz9+ by} g+ 6yh49+ 4449+ 4409), . (8) 
N an(2)=6N (23+ 215 + 2afy + Zagat 2st 215 + Gye 4 ByYbyg+ 2yhq+ 24, 
+ Gib + Gif gt Spy Spy9+ Sprott 4yoo) - - » « . (9) 


The logarithm of a general term of the configurational partition function is 


5N 
=> {wlnw+(1—a) In(1—x)}— 5 Zvah Ing, — 2 (ON ap(1) +N 43(2)}. 
(10) 


The seventeen probabilities 3, ..., %19, %12,--++, %o9 are chosen as 
independent variables and the remaining three ¢,, %,, and yb, are given by 
(5), (6) and (7) respectively. Introduce the symbols 7, €, «, 8, 6 defined by 
exp (—w,/kT')=n (w,>0), exp (w,/kT)=n (wy<0); exp (—w/kT)=E ; 
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pi =a8 2 b2—=P* ; an/B=4. (« and f are real and positive.) Now differ- 
entiate (10) with respect to the independent variables 


Pro= 7088S 

p= F40°B Pyg= Sta BN, 

py=F40H% Py F405 BP 

b= §8a4 6? Ps S9a2B yh 

pe = Sat B? PKS CSL Oe 9 eee eset OI) 

f= E80. 84 P= ESat BA 

pye= F'n? ft Pig = Feat BA 

Po= F833 Py o=°E8a3 BF, 

Po = F803 3 Poo= 195803 BF 4 
To obtain the conjugate phases in the disordered case a differentiation 
with respect to w is necessary, 


Py =a8(1—x)5/ax5, . 7 . a ‘3 : 4 (12) 


This also gives the equilibrium value of s in the ordered case and the 
relative magnetization in the Ising Model. Substituting (11) in (5), (6) 
and (7) and eliminating 8 and 4%,,, 


afl—ax)=frofs(Aifs) + - - . ~ (18) 
where 
Fl P= + 5yPEAh+ (2AE8-+ 8429) $?+ (Gy°Eo+ 4yE%) G5 
+(n°E8+ 4768) $4 + E49", Mee eae eS 3 .. ey (14) 
SAP)=P A$”); Pee: (15) 


SP) =U + 6th + (Bn AEP+ 129268) hb? + (L2FES + 87E9) $? 


DA (8777s2a 127°) p4-+ 6netd?® 4%, peewee. are 1G) 
Ss(¢)=$f3(¢ 7). SB ph sete RS es arcane oe YD 
Using (12), z/(1—zx) can be eliminated to give 

ANOS SSS PCE atte Ree etary Sd @E:)) 


(18) is an equation in the single unknown ¢ and gives ¢ for any pair of 
values of 7 and €; the conjugate phases are then given by (13). 

Choosing the value 0-5 for p (18) was solved numerically for three 
values of 7, 0-77, 0:79 and 0-81, and the corresponding values of x were then 
calculated from (13). 

The critical temperature for p=0-5 is given by 


Ne=exp (—w,/kT,.)=0°8344, 
which was calculated from the equation given in the next section. In 
table 1 the ratios 7/7, are compared with those obtained from the present 


quasi-chemical approximation with p=0 (when it becomes identical with 
the quasi-chemical ‘ pair’ approximation) and the zeroth-order approxi- 


mation. 
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The zeroth-order method for the disordered case becomes accurate when 
the interaction energies are independent of the distance between the 
systems and it is thus to be expected that the effect of fluctuations will be 
diminished by a reinforcing second interaction. This is confirmed by 
table 1 which shows that the 2a—1, 7/7, curve for p=0-5 lies roughly 
half-way between that for p=0 and the zeroth-order curve. 


Table 1 


‘ 


T/T, (quasi-chemical) 
2a—1 (=s) |-#-—_—_——__ 1/T,, (zeroth-order) 
p=0-°5 


0-866 : 0-693 0-657 
0-791 0-768 0-736 
0-654 0-880 0-859 0-836 


§3. THE CriTicAL TEMPERATURE 


At the critical temperature the conjugate phases coalesce at x=4, 
so that, from (13) d=1 and (18) becomes an identity. The critical 
temperature can be obtained by putting ¢=1-++6 in (18) and equating the 
coefficient of 5 to zero. The resulting polynomial equation in 7 and € is 

£16 (67)4 +4 4875+ 21 678+ 120y7-+ 2478) 
+ €14(4874-+ 224n°-+ 2887'+ 5607+ 19278) 
+ €12(— 24n*+ 447° + 9678+ 2727+ 5528+ 108°) 
+-€19(— 80n4—3367>— 1927+ 11278-+ 4329) 
+ £8( — 2n*?— 72y3?— 32n*— 32478 — 2078+ 16671°) 
+ €6(—32n?— 483— 128y4— 8078+ 641°) 


+ &4(—4n—100n>°— 88n’+ 567n")+ 6y!2—24y8=0. . (19) 
Where €=1 this reduces to 
2n(3y—2)(1-7n =O 2 ee (20) 


The only positive solution of (20) is 72/3 whichis, as would be expected, 
the pair approximation result for z=6. 

When the first interaction is zero the correlation between sub-lattices 
a and b is removed and the simple cubic lattice decomposes into a 
pair of face-centred cubic lattices. If the method discussed here is to be 
consistent the critical temperature must, as 7 tends to one, approach a 


value close to the pair approximation critical value for z=12. Putting 
n==1 (19) becomes 


(4628+ 64¢° 28¢418)(98+ 16264 6€4+1)=0. . . . (21) 
The first bracket on the Jeft-hand side gives the only positive root and on 
equating it to zero the critical equation for the face-centred lattice using 


a group with n=6, m=12 is obtained. This yields a value 0-8300 for es 
which is quite near the pair approximation’s 10/12 or 0-8333. 
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For purposes of comparison the critical equation for the plane square 
lattice in which z,=z,—4 has been deduced. A group of a central site 
with its four nearest neighbours, giving M,=m,—=4 and twelve basic 
configurations, was used and the following critical equation obtained : 


£°(2m9) + E9(LOn*+ 813+ 2m?) + €4( 42+ 493+ 2074 + 2075) 
+E —1—4n—9n°— 54+ 4n°+ 1198)—4y 34 2yt=0. . . (22) 
Where €=1 this reduces to 


een let y)) == One) Gage) oes = (23) 
whose only positive root is 4, confirming that for w,=0 the ordinary 
quasi-chemical pair approximation is regained. 

The two equivalent sub-lattices are also of the plane square type so that 


the accurate critical values of 7 where £=1 and € where 7=1 are both 
0-414 (Kramers and Wannier 1941). Where 7=1 (22) becomes 

Mee 467 — 1)(S 610. 2 a. (24) 
The only positive root is given by the first bracket on the left-hand side 
and on equating it to zero the critical equation for the ‘ quadruplet 
approximation ’ for z=4 is obtained. Where w,=0 the critical value of 
€ is thus 0-486 while that of 7 with w,=0 is 0-500 so that the effects of 
first and second interactions appear reasonably well-balanced in this 
approximation. 


—§4. CrivicAL TEMPERATURE AND ENERGY VALUES 


Pairs of values of 7 and € satisfying (19) and (22) were obtained by 
giving 7 a definite value and finding the positive root of the resulting 
polynomial equation in €. Since there was never more than one change of 
sign in the coefficients of this equation the positive root obtained was 
unique by Descartes’ Rule of Signs. Values of In 2/In7,~1 are plotted 
against p (=In é,/Iny,) in fig. 1. Where w,>0 this gives a curve of 
(kT, In 2)/w, against w,/w,. 

If w, is kept constant the critical temperature, 7',, in both the simple 
cubic and the plane square cases decreases smoothly to zero as w, decreases. 
The curves thus resemble that derived by Domb and Potts (1952) for the 
plane square lattice, using a series approximation, and show no sign of the 
‘kinking ’ observable in that derived by Chang for the simple cubic lattice 
(1937). Chang used a quasi-chemical approximation different from that 
employed here. 

The curve for the simple cubic lattice meets the p axis at —I/4; as 
7’, approaches zero the critical values of 4; and 15 approach those of #, 
and %,,. This corresponds to the fact that where p<—1/4 an ordered 
state with A and B systems in alternate layers is more stable at the 
absolute zero than pure A or B phases. 

Where w,<0, p=—vw,/w, and the ordinate in fig. 1 represents 
(kT', In 2)/—w,, T, being now the critical temperature for long-range order 
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with respect to the face-centred cubic sub-lattices @ and 6. The 
ordered state at 70 when p<—+ now consists of rows of A and B 
arranged so that each row of A is surrounded by four rows of B and 
vice versa; this was shown by Luttinger (1951) to be the state of least 
energy where —w,>—}w,>0. The fact that the 7’, curve cuts the axis 
at the correct value of p represents a further check on the consistency of 
the approximation used here. 

The 7’, curve for the plane square lattice cuts the axis at —}. Where 
p<—#anarrangement of alternate rows of A and B is the most stable state 
at absolute zero for both positive and negative values of w,. 


PLANE SQUARE SIMPLE CUBIC 
LATTICE LATTICE 


— 0°25 0 = ; : 
=p 0-25 0-5 


Critical temperature curves. 


For p>0 the significance of the critical temperature results is best seen 
by plotting 7,,°/7',, T° being the zeroth order critical temperature, against 
%2p/(z,+22p). The quantity zyp/(z;+22p) has the value 0 where w,=0, 
1 where w,=0 and } where z,w,=z,W}. 

kT,—=w,|In ne and T',°=}w,(2,+2sp)/k so that for the simple cubic 
lattice 

PoP = 3(ln yt 2 Ine) ee) 
The critical values of N4,(1) and V,;(2) may be obtained from (8), (9) 
and (11). Defining W as the critical value of the configurational energy 


W/RkT = {NV 43 (1) In oe +N ap(2) In nad (EN tes (26) 

The zeroth-order critical energy W® is given by 
, W°=N(z,w,+29W,)/4. ne ctainde etch Mate ea aah 
Since W9/ EPs LB iar neue PUT Pog crue. eee 


(WRT) (Wi ECT Vex abe? Wamaee an oe ero 
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Curves of 7,°/T, and RT/2W are plotted against 2p/(1+2p) in fig. 2 
and both show that the effect of fluctuations is diminished when there & 
reinforcing first and second interactions. It appears from these results 
and from those given in § 2 that the differences between values calculated 
for appreciable p and those calculated for p=0 are a considerable fraction 
of the differences between the latter and the zeroth order values e.g. for 
p=1/8, corresponding to an inverse sixth power law, (W°/RT 0/( W/RT’) 
has diminished from 1-027, its value at p= 0,5to 1-016. ; : 


Fig. 2 
0-3 
0-2 
0-4 
re) 
0-2 0-4 0-6 0-8 {-O 
2p 
I+2e 
Critical temperature and energy curves. 
TO Ww? |W 
1) ~—-1l. 2 —— /—— —l\. 
() 7-1 (2) 10 {ere Fer i} 


W® is the total configurational energy increase from 7=0 to T’'= so 
that W/W? is the fraction of this increase which has taken place at T=. 
W/W? is plotted in fig. 3 against 2p/(1+-2p) over the whole range from —1 
to 1 in which transitions of the type defined by (19) are possible. The 
values for p>0 again illustrate that the division of the configurational 
energy between reinforcing first and second neighbour interactions 
diminishes fluctuation effects. The opposite occurs where p<0 and 
W/W® decreases rapidly to zero at 2p/(1+-2p)=—1 or p=—I1/4. Van 
Santen (1951) considered that low values of W/W® for appreciable further 
neighbour interaction are necessary to account for the observed magnitudes 
of latent heat of melting for ionic solids. 

Chang (1941) showed, by means of Kirkwood’s method, that an increase 
in critical short-range order must occur for negative values of p. Defining 
the short-range order 

__ Naa()+Nnn(l)—Nan(t) 


a . (30) 


74 G. M. Bell on the Second 


it is found from (19) that, at the critical point for the simple cubic lattice, 
where 

p=0 ,  o=0-2000, 

p=—0-0745, 00-2330, 

p=—0:1421,  o=0-2811. 


zl O 
20 
I+2p9 


Critical energy ratios. 


§5. OTHER QuUASI-CHEMICAL METHODS 

Chang (1937) applied the Bethe method to the simple cubic lattice for 
the p<0 case using a basic group consisting of a central site, an inner layer 
of six first neighbours and an outer layer of their eighteen first neighbours. 
Hence, n=25, m,=36 and m,=72. A complete consideration of this 
group would give the next approximation to that employed above but 
would be extremely laborious. However, Chang only considered 24 of 
the 72 second neighbour links giving an effective z, of 4 instead of 12. 

The ‘ quadruplet ’’ approximation with n=m,=4, m,=—2 has been 
applied to the simple and body-centred cubic lattices for the positive 
values of p corresponding to an inverse sixth power law by Guggenheim 
and McGlashan (1951 b). 

These authors recognize that, since the effective value of z, is always 4z,, 


the approximation is inconsistent when z,4}2z,. The critical equation is, 
in the notation used here, 


1-2? 48 /e, ae") — 3 fz 0. 8... BTS 
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Where €=1 the quadruplet approximation for first neighbours is regained 
but the results of the inconsistency appear when y is put equal to 1 and 
(31) becomes 
| Pde —laai=0, . 2... (32) 
For z,=4 the effective value of z, is 2 and there is no positive solution 
of (32). For z,=6 the effective z, is 3 and (32) gives the solution fa) 3 
For z,=8 the effective z, is 4 and (32) gives £=1/2. Now the body- 
centred cubic lattice can be divided into two equivalent simple cubic 
sub-lattices, such that any site and its second neighbours lie in the same 
sub-lattice. The limiting value of é as w, tends to zero ought, therefore, 
to be near 2/3. For all values of z, the critical temperature given by 
(31) will go to zero at p=—1, though the correct values are p=—1/4 and 
p=—2/3 for the simple and body-centred cases respectively (Luttinger 
1951). It is thus obvious that the quadruplet method, with its low 
effective z,, underestimates the second neighbour effects and, in fact, 
Guggenheim and McGlashan come to the conclusion that these are of 
trivial importance in relation to measurable quantities. However, 
differences between results calculated on the quadruplet approximation 
with and without the second neighbour perturbation do not always seem 
to be negligible in comparison with the differences between the latter and 
the zeroth-order results (Guggenheim and McGlashan 1951 b, table 3). 
Hill (1950) uses the quadruplet approximation for the plane square 
lattice and to correct for the inconsistency he regards the factor & as 


exp (—2z,w,/z,kT), 


_ instead of exp (—w,/kT'). However, in the limit w,=—0 there is again no 
finite critical temperature and, although Hill does not suggest using the 
approximation in the region of large p, this fact does not indicate that the 
correction is very effective in the plane square case. It is of some interest 
to apply our criteria to the modified quadruplet approximation for the 
simple and body-centred cubic cases. For 2,=6, z.=12, w,=0, 
exp (—w,/kT’,)=0-7598, a value which corresponds to an effective 2, of 
about 8; 7’, becomes zero at the correct value, —1/4, ofp. Forz,=8, zg 6, 
w,=0, exp (—w,/kT’.)=0-6300, a value which corresponds to an effective 
z. of between 5 and 6; 7’, again becomes zero at the correct value of p, 
which is now —2/3. The correction thus results in a considerable improve- 
ment in these two cases. 

The quadruplet approximation has not been applied to the face-centred 
cubic lattice where, since z,=6—}z,, it gives the right ratio of second to 
first neighbour pairs. This lattice is composed of four simple cubic sub- 
lattices such that a given site and its second nearest neighbours all lie in the 
game sub-lattice. Thus, as w, tends to zero the critical temperature should 
tend to that of the simple cubic lattice. Putting 2; =12 in (32) we obtain 
the positive root 2/3. 

Where w,>0, w,<0 and p<—lI pure A and B phases (ferromagnetic 
state) are less stable at 70 than a state in which there is a completely 
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ordered arrangement of A and B on each sub-lattice, the sub-lattices 
being thus uncorrelated (Smart 1952, table 2). 

Equation (31), in making the 7’, curve cut the p axis at —1, is thus 
correct for the w,>0 case on the face- centred cubic lattice. Hence, by 
this criterion ana by that of the limiting values at w,=0 and w,=0, the 
quadruplet approximation is a consistent one for the disordered phase 
separation or ferromagnetic case on the face-centred lattice. Since this 
lattice cannot be decomposed into two equivalent sub-lattices there is no 
simple connection between the equilibrium states for w,>0 and those for 
w, <0, these latter being rather complicated (Anderson 1950, Smart 1952). 
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ABSTRACT 


A coincidence experiment is described in which an attempt is made to 
detect fissions of uranium nuclei resulting from the capture of negative 
#-mesons in the cosmic radiation. The result shows that fission is not 
the most probable process following meson capture, and an upper limit 
of 0-25 is set for the probability of this fission process occurring. 
A interpretation of the result is suggested in terms of the competition 
between fission and neutron emission in the protoactinium isotope 2°8Pa 
which is formed in the capture process. 


§1. [yTRODUCTION 


THE processes which occur when a negative u-meson is captured by a 
nucleus are usually explained in terms of the so-called charge-exchange 
reaction, in. which the meson interacts with one of the protons of the 
nucleus, giving rise to a neutron and a neutral particle which may be a 
neutrino 


p-+P=N-+». 


The amount of excitation energy imparted to the resultant nucleus, in a 
charge-exchange reaction in which a neutrino is emitted, has been 
calculated by Tiomno and Wheéler (1949) using three different nuclear 
models. The three calculations give results of the same order of 
magnitude and the free particle model, which is most applicable in the 
present case, predicts a most probable excitation of ~ 14 Mev and a 
maximum of ~ 22 Mev. Rosenbluth (1949) finds a value which is almost 
exactly in agreement with Tiomno and Wheeler, but Taketani, Nakamura 
and Sasaki (1949) calculate the nuclear excitation to be about 6 Mev 
for large Z and about 14 mev for small Z. Lang (1952) has given reasons 
for concluding that the excitation energy is about 1-5 times that predicted 
by Tiomno and Wheeler. 

The excitation of the nucleus to the predicted extent would be expected 
to result in the emission of neutrons, and these have been observed, from 
lead, by Sard et al. (1951), Groetzinger and McClure (1950) and Sard 
and Crouch (1952). The number of neutrons emitted is about two, which 
is in reasonable agreement with theoretical prediction based upon the 


* Communicated by the Authors. 
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evaporation formula of Weisskopf. George and Evans (1951), using 
photographic emulsions underground, have shown that only about 9% 
of the captured mesons produced disintegrations accompanied by the 
emission of charged particles, and that 90% of these visible stars had one 
prong and 10% two prongs. This, according to the calculations of 
Fujimoto and Yamaguchi, is consistent with a nuclear excitation of about 
15 mev. The charge-exchange reaction as given above is therefore 
supported by a considerable body of experimental evidence. 

Wheeler (1949) has pointed out that the capture of a ~~-meson by a 
uranium nucleus might be expected to lead to fission by one of two 
possible processes: (a) The meson is first trapped and moves in stable 
orbits round the nucleus. The calculated energy released in the Is—2s 
transition is 7 Mev uncorrected and 5-6 Mev if allowance is made for the 
dumb-bell shape of the nucleus immediately prior to fission. If this 
energy, instead of being radiated as a hard x-ray, could be transferred 
to the nucleus, fission might result. (b) In heavy elements, the meson 
is captured by the nucleus before it decays. As has been shown 
above the excitation energy imparted to the nucleus is of the right 
magnitude to induce fission. 

The observation of fission induced by the capture of «~-mesons, if it 
occurs, would be intrinsically difficult, since it would have to be observed 
in the presence of the natural «-activity and also of spontaneous fission. 
The rate of capture by uranium of »~-mesons from the cosmic radiation 
is calculated to be about 7 x 10~? per gramme per hour, while the measured 
spontaneous fission rate is 24:8+0-9 per gramme per hour, and the 
a-emission is about 105 per gramme per hour. Simple calculations 
showed, however, that coincidences between the passage of mesons 
through counter trays and the occurrence of fissions in an ionization 
chamber could be observed, and that the background of random coinci- 
dences would be negligibly small. If the majority of the captured 
mesons had given rise to fissions, the number observed, though small, 
would have been statistically significant. In fact, the number of events 
recorded was very much below the calculated maximum, which indicates 
that fission is not the most probable consequence of u~-meson capture 
by a uranium nucleus. 


§2. EXPERIMENTAL ARRANGEMENT 


A fuller account of the apparatus and of the precautions taken to ensure 
reliability has been given elsewhere (Galbraith and Whitehouse 1952) 
and only the essentials will be described here. 

The general arrangement of the trays of Geiger-Miiller counters is 
shown in fig. 1. The upper trays A and B, which had areas of 4900 em? 
and 3060 cm? respectively, detected the passage of penetrating charged 
particles into the apparatus. The tray D was connected in anti-coinci- 
dence, and, although it could not be made as efficient as is usual, owing 
to the necessity of supports for the ionization chamber, served to reduce 
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the background of random coincidences by a large factor. The flat, 
cylindrical ionization chamber C contained about 8 grammes of uranium 
in the form of thin foils of U;0,. The measured rate of spontaneous 
fission in this chamber was 115 per hour, and this, in conjunction with 
the known spontaneous fission rate ($1), gave the effective weight of 
uranium, from which fissions could be observed, as 4-7+-0-5 grammes. 
The events of direct interest are the passages of charged particles through 
A and B but not through D, accompanied by a fission in coincidence. 
Such events are denoted by (AB-D : (). 

The function of the lead between the two upper trays A and B was to 
prevent the production of false (AB—D : C) events due to the discharge 
of the counters by y-rays emitted in spontaneous fission. The thickness of 
lead was chosen to decrease the y-counting rate in A by so large a factor 
that this effect was excluded. Subsidiary experiments, in which a 
Po-«—Be neutron source was used to increase the fission rate in the 
chamber, were carried out to verify this point. 


Fig. 1 


Arrangement of counter trays and ionization chamber. A, B and D indicate 
the central planes of the trays of counters. C indicates the mean position 
of the uranium foils. S indicates the approximate position of the shower 
tray. 


The tray S marked in fig. 1, which consisted of eight counters coupled 
alternately in two groups of four, was used in part of the experiment to 
ensure that the apparatus was not sensitive to dense electron showers. 

A block diagram of the electronic apparatus is given in fig. 2. When 
an (AB: C) event occurred the fission pulse was displayed on a cathode 
ray oscilloscope and photographed, as a precaution against spurious 
events. The other events (AB-D:C) were recorded on pen recorders. 
The resolving times in the coincidence unit were set at one micro-second, 
and subsidiary experiments were carried out to ensure that the pulses 
from all the units, particularly those from the ionization chamber, arrived 
within this time, and that no coincidences were lost. Daily checks were 
made of the reliability of the whole equipment. The spontaneous fission 
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rate and the (AB-D) rate were recorded continuously as a further check 
on the operation of the equipment. The whole apparatus was driven from 
a motor generator to reduce the possibility of spurious counts due to 
surges in the mains. a3 

The apparatus was first set up at ground level under a jj in. steel 
roof, and later moved below a concrete roof, 3 ft. thick. 


Fig. 2 
OSCILLOSCOPE OSC. TRIGGER PULSE 
CLOCK AND 
CAMERA C PULSE 
‘A TRAY AL ABC 


NORMALLY USED 
To RECORD ABC. 


AMP. 


a ri 
D 


MAIN | SELECT 


TEST = 
Purse, D TRAY 
ie ABP [SCALEARER 


mag, ABC: SHOWER 
POOR) coNcipENCES 


— | COINC. | SHOWER 
SHOWER TRAY pcs UNITS re 


(USED IN PART OF EXPERIMENT ONLY) 


Block diagram of electronic apparatus (P.O.R.=Post Office Register ; 
P.R.=pen recorder. The ‘ select output’ scaler and recorder, normally 
used to record ABC coincidences, could be used to measure the A, B, D, 
AB, ABD and AB-D rates). 


§ 3. CALCULATED RATES 
3.1. Rate of Capture of u--mesons 

The flux of penetrating charged particles at sea level was not measured. 
in this experiment. Instead the result of Greisen (1942, as reported by 
Rossi 1948), who found a value of 46 particles per cm? per hour, was 
used. Of these about 99% are u-mesons and 1°% fast protons (Mylroi 
and Wilson 1951). Further, 45% of the mesons are negatively charged 
(Owen and Wilson 1951). A calculation of the geometrical efficiency of 
our counter trays shows that 90% of the particles which pass through 
the uranium also pass through the trays. Allowing for the combined 
detection efficiency of the trays A and B (93%), we estimate that the 
effective flux of .~-mesons through the uranium is 17 per cm? per hour. 

The fraction of 4.-mesons stopped by uranium may be calculated from 
the absorption in lead, which has been measured several times (Koenig 
1946, Rossi 1948, Kraushaar 1949, York 1952). We have used the most 
conservative value, that of Koenig, who found that 1 g per cm? of lead 
stopped 0-033°% of the mesons. Assuming the mass absorption of the 
two elements to be the same, and the effective mass of uranium to be 
4:7 g (§ 2) the rate of absorption of ~-mesons, and hence the maximum 
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possible rate of u-meson fission, is 27 per 1000 hours. For that part of 
the experiment which was carried out under concrete the meson flux was 
attenuated by about 15%, but the energy spectrum was very little changed, 
and the proportional absorption of mesons was the same as above ground. 


3.2. Proton Reactions 

The flux of fast protons with momentum great enough to penetrate 
the lead between the trays and enter the chamber (~ 10° ev) is about 
0-2 per cm? per hour (Mylroi and Wilson loc. cit.). If we assume that 
the cross-section for interaction is the geometrical cross-section of the 
nucleus and that every interaction results in a fission, we obtain, for the 
_ present apparatus at ground level, a maximum rate of 7 proton-induced 
fissions per 1000 hours. Below the concrete roof the proton flux is 
reduced by a factor of about 5, if an attenuation length of 140 g per cm? 
is assumed (Mylroi and Wilson Joc. cit.), and the number of proton-induced 
fissions is therefore hardly observable. 


3.3. Other Mechanisms 


There are a number of indirect processes by which fast protons may 
initiate (AB—D : C) or (ABD: C) events. Most of these depend upon the 
production of a star in the material of the apparatus and the subsequent 
induction of a fission by a secondary particle from the star. Approximate 
calculations indicate that the fission rates from these processes are much 
smaller than those from u-meson capture above ground, while below the 
concrete they are virtually absent. 

Another type of reaction which might conceivably lead to fission is 
the direct electromagnetic interaction which occurs when a fast charged 
particle passes close to the nucleus. The results of George and Evans 
(1950) indicate that the cross-section for this interaction is ~ 10~?° em? 
per nucleon, and a simple calculation shows that fissions caused by this 
effect would be too infrequent to be detected by the present apparatus. 


3.4. Accidental Rates 
The accidental rates can be calculated from the usual formula 
(R,=—27R,R,), in which R, is the accidental rate, 7 is the resolving time 
of the coincidence unit (1 sec), R, is the spontaneous fission rate in 
the chamber and R, is either the (ABD) or the (AB—-D) rate. The 
calculated accidental rates are shown in the table. | 


.§4. RESULTS 


The results of all the runs are recorded in the table. The following 
conclusions may be drawn from them. 

RunsIandII. Together these indicate that the rate of meson-induced 
fission is significantly below the calculated rate, but the interpretation is 
complicated by the fact that the number of (ABD: C) events is above the 
calculated random rate. This may be due to the indirect processes 
mentioned in § 3.3. 
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Run II. No coincidences were observed between (AB: C) events and 
electron showers. 

Runs III and IV. These are satisfactory control runs, in which the 
rates are in agreement with the calculated random rates. The count of 
4(ABD : C) events in run IV is not large enough to justify any definite 
conclusions about the fission of lead by cosmic ray particles, but it would 
not be inconsistent with fission by fast protons at a rate determined by 
the geometrical cross-section. 

Run V. The number of (AB—D:C) events is so much less than the 
calculated rate of absorption of .~-mesons that the experiment provides 
no evidence for fission by z-mesons in uranium. The statistical accuracy 
is such that the upper limit of the probability of fission following the 
capture of a .~-meson in *°8U can be set at 0-25. 


§5. Discussion 


At first sight the results of this experiment seem to be at variance with 
predictions based upon the charge-exchange reaction. It seems probable, 
however, that an explanation for the result can be made in the following 
way. 

Experimental evidence has accumulated which indicates that fission 
competes unfavourably with neutron emission at moderate excitation 
energies. It has been shown that the photofission cross-section of ?88U 
is only a fraction of the cross-section for neutron emission at energies of 
the order of 20 Mev (Charbonnier, Scherrer and WaAffler 1949, Goward, 
Jones, Watson and Lees 1951) in spite of the fact that the photofission 
threshold, 5-1 Mev, is lower than the threshold for the (y, n) reaction, 
5-9 mev (Koch, McElhinney and Gasteiger 1950, Anderson and Duffield 
1952). Thus it appears probable that neutron emission is favoured with 
respect to fission even when the fission barrier is less than the binding 
energy of the neutron. 

Price and Kerst (1952) estimate, from their observations of the photo- 
neutron yield, that only about half of the °°U nuclei undergo fission when 
bombarded with 22 Mey y-rays, the rest losing most of the excitation 
energy by neutron emission. For *Th they find the probability of 
fission occurring is only 0-24 for the same excitation energy. Levenger 
and Bethe (1952) show that the photo-neutron yield for **U observed 
experimentally by other workers is consistent with a value of 0-23 for 
the probability that fission occurs in place of neutron emission. We would 
expect that the probability of fission occurring in the compound nucleus 
288Pa, which is formed when a p-meson is captured by **°U, will be even 
smaller, since its Z?/A value is less than that for °°U. The upper limit 
of 0:25 observed in the present experiment for the probability of fission 
occurring would then appear to be in reasonable agreement with the 
above considerations. 

If this explanation of the failure to observe j.-meson fission is correct, 
there is no conflict between the result of the present experiment and other 


work on p-meson capture. 
G2 
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SUMMARY 


Nuclear-research emulsion has been exposed to the cosmic radiation 
under various thicknesses of ice. Three events interpreted as the decay 
of r-mesons have been observed, and measurements on their tracks show 
that these particles decay into three 7-mesons of total kinetic energy 
(73-+6) Mev and have mass (965+12)m,. They are found with a frequency 
of about 0-0008 z-mesons per star in emulsion exposed under ice. The 
emulsion was examined for nuclear disintegrations produced by heavy 
mesons, and for z7-mesons decaying in flight, but no such events were 
found. The life time is shown to be longer than 10-9 sec. 


§1. INTRODUCTION 


In 1949, Brown, Camerini, Muirhead, Powell and Ritson (1949) observed 
a single event which indicated the existence of a particle having a mass 
about a thousand times that of the electron, and decaying into three 
singly-charged particles. This particle was called the 7-meson, and its 
existence was confirmed when, during an investigation of the attenuation 
of cosmic radiation in ice, Harding (1950) found two events interpreted 
as the decay of 7-mesons (numbers 1 and 2 in this paper). 

A second batch of plates has been exposed in the ice and examined for 
further evidence of t-mesons. One more (number 3) was found (Hodgson 
1951). 

The measurements on all the 7-mesons so far observed at Imperial 
College have been repeated and extended, and the results are given in 
this paper together with the evidence about the mass, mode of decay, 
frequency of production, and life time of these particles. The plates 
have also been examined for evidence of slow heavy mesons producing 
nuclear disintegrations. 


§2. EXPERIMENTAL PROCEDURE 


Boxes of plates coated with Ilford G5 nuclear-research emulsion were 
exposed vertically in a steel pipe sunk into the ice on the J ungfraujoch 
plateau (see Harding 1951). Particles 1 and 2 were found in 200-y-thick 
emulsion, number 3 in 400-~ emulsion. 


* Communicated by E. P. George. 
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§3. THE Mass oF THE t-MESON 


(a) Gap measurements on the tracks of the t-mesons 

The results of gap measurements on the tracks of particles 1 and 2 
have been published (Hodgson 1950). However, as the tracks had faded 
considerably, it was not possible to find an accurate value of the mass 
of the r-meson from them. The track of the remaining t-meson dips 
too steeply for any measurements to be made by this method. 


(b) Multiple scattering of the t~mesons 

Measurement of the multiple scattering on the track of particle 2 gave 
a value of (1870-+500)m, for its mass. This is rather high, but it is not 
inconsistent with the more accurate value found below. The tracks of 
the other 7-mesons were not suitable for scattering measurements. 


(c) The energy balance of the decay process 

The most accurate value of the mass of the r-meson is obtained by 
assuming the decay scheme and adding the rest masses of the secondary 
particles to the mass equivalent of their kinetic energies (see §4.) The 
best value found in this way is m,=(965+12)m,, in excellent agreement 
with the figure of (966--8)m, given by Fowler, Menon, Powell and Rochat 
(1951). In both these investigations, the mass of the secondary 7-mesons 
was assumed to be (274+2)m, (Alvarez 1950). Since this paper was 
written another 7-meson has been reported by Ceccarelli, Dallaporta, 
Merlin and Rostagni (1952) who find a mass of (991+15)m, which is 
consistent with our value. 


§4. Tor Drcay OF THE t-MESON 


In all the examples of 7t-meson decay observed so far, the tracks of 
the three secondary particles are coplanar to within 2°, indicating that 
not more than three particles are emitted. The possibility of the 
emission of a neutral fourth particle of low energy or small mass cannot 
be excluded definitely, but it will be disregarded in the analysis to follow. 

The grain densities of all the tracks of secondary particles suitable 
for measurement were found, and they are given in table 1 as a multiple 
of the minimum grain density of singly-charged particles. Owing to 
fading, the values quoted for the first two events are not as reliable as 
those for the third. 

The relative masses of the secondary particles were found from their 
directions of emission and the grain densities of their tracks by the 
methods of Fowler et al. (1951), and they are given in table 1. These 
results show that the secondary particles are of equal mass to an accuracy 
of about 10 per cent. 

The absolute values of the masses of those secondaries having a 
sufficient length of track in the emulsion were determined from measure- 
ments of grain density and multiple scattering, and the results are shown 
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in table 1. These masses are in the range 200 to 400m,, indicating that 
the secondary particles are - or a-mesons. But since one of the 
secondaries in the event observed by Brown et al. (1949) was identified 
as a m-meson, and since it has been shown above that the secondary 
particles are all of the same mass, we can conclude that all the secondary 
particles are 7-mesons. Provided that all the 7-mesons observed have 
the same decay scheme, the probability that one or two of the secondaries 
are z-mesons is negligible. 


Table 1 
7-Meson 
1 4 3 
I/Imin of Secondaries 1-62 -+0-30 — 1-90 40:25 
2-90 -+0-90 2:36 -+0-20 2:00-40-15 
3°52 -+0-50 —— 4-72+0-70 
Relative Masses of 1-00 — 1-00 
Secondaries 0:96 0-25 — 0-93 0-12 
1-05 -L0-20 — 1:09 +0-17 
Absolute Masses of — a 260-++100 
Secondaries (m,) — — 330 +140 
300 +70 a 260-110 
Relative Momenta of 1-00 = 1-00 
Secondaries* 0:-59-+0-18 — 0:95 -+0-10 
0:-51+0-11 — 0-49 +0-06 
Relative Momenta of 1-00 1-00 1-00 
Secondariest 0-57 0-96 0-88 
0-53 0-91 0:54 
Energies of Secondaries 45-0+12 25 +5 33-5 +5f 
in Mev ieee) 2414 Ss) Osed 
13-3+ 3 20 +4 9-0-+1-5 
Total Kinetic Energy 76+15 69-+8 73-547 


of Secondaries in Mev 


* From grain densities. 


+ From angles of emission. 
+M. G. K. Menon of Bristol University has measured the energies of the 
secondary particles using the scattering technique, and he finds them to be 


33-+10, 27-6 and 8-1 mev respectively. 


Knowing the secondary particles to be 7-mesons, their energies were 
found from the grain densities of their tracks ; they are given in table l. 
The relative momenta of the secondaries are also given, and they agree 
well with the values calculated from the angles of emission. 


88 A. J. Herz, P. E. Hodgson and R. M. Tennent on the 


The values for the total kinetic energy of the secondary particles agree 
within the limits of experimental error, supporting the view that 7-mesons 
decay into just three particles. The mean value of the total kinetic 
energy of the secondaries is (73.6) Mev, in agreement with that given 
by Fowler et al. (1951). 

In the case of the second t-meson, two of whose secondaries were too 
short even for grain-density measurements, the total kinetic energy of 
the secondary particles was found from the angles of emission and the 
energy of the third particle. 


§5. THe Rats oF PRODUCTION OF 7-MESONS 


During his investigation of the attenuation of the cosmic radiation in 
ice, Harding (1950, 1951) observed two 7-mesons among 3673 stars of 
three or more tracks in plates exposed in the ice at depths between 
0 and 9 metres, which gives a ratio of 7-mesons to stars of 0-000579 S001. 
In the second exposure, made at depths between 2-5 and 3-5 metres of 
ice, one t-meson was found among 1300 stars and 170 7—y decays and 
meson-induced stars, giving r-meson frequencies of 0-00087 950g per star 
and 0-006+0°14 ner meson of the type given above. 

Fowler et al. (1951) examined plates exposed under 30 cm of lead, and 
found one 7-meson accompanied by 750 z-mesons, that is 0-0013%(°098 
per 7-meson. 

The errors given here were computed from the fiducial limits given by 
Regener (1951), and it can be seen that the frequencies under ice and 
under lead are of the same order of magnitude. However, the difference 
between the figures suggests that 7-mesons are more frequently found 
under ice than under lead, and if this were confirmed by further 
experimental evidence it would provide support for the hypothesis 
(Hodgson 1951) that 7-mesons are more likely to be produced in hydrogen 
than in other elements. This may be understood if it is assumed that 
T-mesons, or their parents if they are secondary particles, if created 
in nucleon-nucleon collisions in larger nuclei, have a high probability 
of reabsorption in the nucleus in which they were produced. 

The 7-mesons made angles of 85°, 169° and 133° respectively with the 
upward vertical when they entered the emulsion. These data are 
insufficient to allow positive conclusions about the zenith-angular 
distribution of the particles to be drawn. 


§6. THe Lirz TIME oF 7-MESONS 
The fact that the frequency with which r-mesons are found to stop in 
plates surrounded by ice is comparable to that in plates surrounded by 
lead, or even greater, is the more surprising because for charged particles 
of short life time the probability of coming to rest increases with the 
density of the medium in which they are slowed down. This probability 
has been calculated for particles of various energies, using the methods 


described by Herz (1951), and an analysis based on the results is given 
in this section. 
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If we assume that all the z7-mesons found came from the ice, we can 
calculate how many had to be present just outside the pipe for every 
one stopping in the plates. At that point the particles would have had 
a kinetic energy of about 55 Mev. 

However, certainly not all the t-mesons were created just outside the 
pipe, and it is therefore of interest to find out how many would have to 
be produced per particle stopped if production took place at some higher 
energy. Now, so far, there has not been any positive experimental 
evidence to suggest that z7-mesons are ejected from nuclear interactions 
with kinetic energies sufficiently low to make them identifiable (Camerini, 
Fowler, Lock and Muirhead 1950, Daniel, Davies, Mulvey and Perkins 
1952). Furthermore, although the accuracy of mass determinations is 
low in cloud-chamber work, it seems safe to say that 7-mesons occur only 
very rarely amongst the non-relativistic decay products of V- and 
K-particles, if at all. 

We shall thus assume for the purpose of this calculation that most of 
the z-mesons are produced with kinetic energies equal to their rest 
energy (about 500 Mev) or greater, and that they are not very abundant 
at that energy (if they were, one would expect a fair number to be produced 
at lower energies also). Table 2 shows how many particles would be 
needed at 500 Mev for every one stopping in the plates; the number 
needed at 55 Mev (just outside the pipe) is given for comparison. The 
figures strongly suggest that the life time is longer than 10~® sec unless 
most of the particles are produced with low energies. 


Table 2 

Attenuation gR ia Life Time in Seconds 
Cross Section Mev 

10-7 10-8 1@-? 10-10 
Negligible 155) 1-003 1:023 1:26 9-76 
Geometrical 55 1-022 1-043 1:28 9-95 e 
Negligible 500 1-068 1-93 713 . 4x Ae 
Geometrical 500 6-68 12-1 4460 PPA See 


A further approach to the problem can be made by comparing the 
experimental upper limit to the number of decays in flight of 7-mesons 
accompanying those stopped with the expected number, assuming various 
life times. Clearly, if more particles are produced than are found to 
stop, the particles which do not survive must disappear either by decay 
in flight or by nuclear interaction. However, no event of either type, 
involving an identifiable r-meson, has yet been found in the eee 
reported here (see also §7). In order to be certaim about this, 0:27 em 
of emulsion was specially examined under high magnification, and, using 
the fiducial limits given by Regener (1951), we can conclude from that 
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alone that the number of decays in flight was not more than 120 per 
z-meson stopped. If we accept the result that no 7-meson decayed in 
flight in the emulsion of the second exposure which was searched specially 
for this type of event, the experimental upper limit becomes 1-9 decays 
in flight per particle stopped. 

Assuming that all 7-mesons are produced with kinetic energy 500 Mey, 
the number expected to decay in flight in the emulsion was calculated. 
The results are given in table 3; they may be in error by a factor of 
about three because of the simplifying assumptions made. Comparison 
of the calculated figures with the experimental upper limit again leads 
to the conclusion that the life time must be longer than 10~° sec. 


Table 3. Expected Number of Decays in Flight of Particles Entering 
the Emulsion with Kinetic Energies between 130 and 490 Mev, per 
7-meson Stopped 


Life Time in sec 10-8 10-9 1g-10 
No nuclear interaction 0:3 250 2-2 x 1029 
Interaction cross section 
taken as geometrical 0-6 900 L3G 


§7. Tue NuciesR INTERACTIONS OF HEAvy MESONS 


Investigations of very-high-energy nuclear interactions by Daniel et al. 
(1952) have shown that few, if any, 7-mesons are produced in them. 
If these particles were produced directly, we should expect them to 
interact strongly with nuclear matter, in which case only those charged 
positively which are repelled by nuclei would be found to decay spon- 
taneously after coming to rest. The negatively charged heavy mesons, 
presumably produced together with the positive ones, will interact with 
nuclei when they come to rest in the emulsion, causing disintegrations. 
It is therefore of great importance to search for more events of the type 
observed by Leprince-Ringuet (1949) and Fry and Lord (1952) which 
can be interpreted as the nuclear interaction of a heavy meson. 

All the stars in the present batch of emulsion were therefore examined, 
but no event was found which could have been produced by a slow 
negative heavy meson. This does not, however, show that such events 
do not occur, for the track of an incoming heavy meson would show little 
multiple scattering and might easily be mistaken for that of a slow 
outgoing proton unless a considerable length of track was available for 
examination. 

If 7-mesons do not interact strongly with nuclear matter, all the events 
could be attributed to the decay of both negative and positive particles. 
This would suggest that 7-mesons are the decay products of still heavier 
particles which are themselves produced in nuclear interactions. 


Production of t-Mesons in Ice sil 
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ABSTRACT 


The image forces on a screw dislocation due to a grain boundary, 
or a free surface on which is a surface film, are evaluated. For the 
case of an oxide layer on aluminium it is shown that there is an 
equilibrium position for the dislocation, at a depth in the metal 
approximately equal to the thickness of the oxide film. The behaviour 
of an edge dislocation is approximately the same. 


§ 1. INTRODUCTION 


Ir has been shown (Eshelby 1951) that there exists a close analogy 
between screw dislocations and electrostatic line charges. In particular, 
the interaction of a screw dislocation with changes in the elastic constants 
of the medium is analogous to the electrostatic problem of a line charge 
in a medium of non-homogeneous dielectric constant. 

The solution of electrostatic problems by the method of images 
(Jeans 1911) can thus be applied directly to the analogous dislocation 
problem. 

§2. GraIn BOUNDARY 


As a model for a grain boundary, consider an infinite isotropic elastic 
medium of shear modulus p,, for x>0, and modulus p, for *<0. Then 
the interaction between a screw dislocation of unit strength at x=a, 
y=0 and the boundary can be represented by the following images. 
For «<0 the stress field is the same as that due to a dislocation of 
strength 2y,/(w e+ ,) at v=a, y=0 in an infinite medium. For «>0 
the stress field is the same as that due to the original unit dislocation 
plus a dislocation of strength (4 3—;)/(jc2-+j41) at the image point x=—a, 
y=0. Thus the dislocation will be attracted to the image (i.e. to the 
boundary) if 2<p,, and repelled if x.>,. This is a long range force 
which varies as a~!. For a real crystal boundary there is an additional 
short range force due to the misfit of the two atomic lattices at the 
boundary (Nabarro 1952). 


§3. Surrace Firm on Free SuRFACE 


Consider a semi-infinite isotropic elastic medium of shear modulus 
#, for «>0 and modulus p, for —b<a<0, with a free surface at x=—b. 


Then a screw dislocation of unit strength at aa, y=0 will induce a 
a eee ee, arene 


* Communicated by F. C. Frank. 
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stress field for x>0 which can be described by an infinite set of images 
at y=0, with values of x, and strengths, as follows :— 


Position of image : Strength of image 
value of a 
—q K 
—(a-+2b) —(1=#2) 
—(a+4b) —(l—k?)« 
—(a+6b) —(1— 12) .c2 
—(a+8b) —(1—?)K3 

etc. ete. 


where k= ((12—1)/(Ha+H1)- 

If wa<p, then the dislocation is attracted to the boundary for all 
values of a. If w,>p, the dislocation is attracted to the boundary 
for large values of a and repelled for small values of a. Thus there will 
be a position where the dislocation is in stable equilibrium, this being 
when the stress due to all the images is zero at the position of the 
dislocation. If this position of equilibrium is ~=d, then, writing n=d/b, 
the condition of zero image stress at v=d is 


n n n n 
0=x—(1—x«?) | at apst es an abiaagt ee lt 


On multiplying this series by «”*! and differentiating a geometric 
progression is obtained. Summing this and then integrating gives as 
the equation for the equilibrium position of the dislocation 


O=K—n(1—e rt | #" (1—t)-1 dt. 
0 


This has been solved numerically for n=0-5, 1, 2 and the corresponding 
values of o/u, are given in the following table :— 


n bolt 
Small 1+2n 
0:5 2-20 
1 3°65 
2 6-76 


For the case of an oxide film on aluminium, this analysis can be 
expected to be approximately true as the elastic anisotropy of aluminium 
is small (Boas 1947). The elastic constants of bulk aluminium oxide 
(Bhimasenachar 1950) are from three to four times the elastic constants 
of aluminium. If the oxide layer on the metal has the same elastic 
constants as the bulk oxide then, by table 2, the equilibrium position 
of a dislocation would be at a depth in the metal approximately equal 
to the thickness of the surface layer. 
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§4. Epae DIsLocAaTION 


An exact analysis (Head, unpublished) of the interaction of an edge 
dislocation with the model grain boundary leads to a rather complicated 
result. However, a good approximation (with a maximum error of 
15%) is that an edge dislocation is influenced by the boundary in 
the same way as a screw dislocation, i.e., as if there were a dislocation 
of strength (u,—p,)/(u2+p) at the image point. The more complicated 
case, of interaction with a surface film, does not appear to be readily 
soluble. The following physical reasoning suggests that there will be 
a position of equilibrium if ~,>,. For when the dislocation is close 
to the surface, the surface film will be effectively infinitely thick and the 
dislocation will be repelled if ».>p,. And if the dislocation is far from 
the surface, then the effect of the surface film will be small and the 
dislocation will be attracted as to a free surface. In between will be 
a position of equilibrium, but its location can only be surmized to be 
approximately that for a screw dislocation. 
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ABSTRACT 


The neutron total cross sections of carbon, aluminium, copper, cadmium 
and lead have been determined as a function of energy in the range 30 to 
139 Mev. The lead total cross section shows a dip at about 60 Mev and 
corresponding dips in the cadmium and copper cross sections appear to 
occur between 30 and 40 Mev. The neutron total cross section of hydrogen 
has been measured as 56-9-+1-8X10-°? cm? at a neutron energy of 
126+2 Mev. 


§1. INTRODUCTION 


NeEvutron total cross sections of several elements at 153 Mev have been 
determined previously (Taylor, Pickavance, Cassels and Randle 1951 a). 
Carbon, aluminium, copper, cadmium and lead have now been studied in 
the energy range 30 to 139 Mev and the cross section of hydrogen has been 
measured at 126+ 2 Mev using the method previously described (Taylor, 
Pickavance, Cassels and Randle 1951 b). Some of the preliminary results 
for lead have already been briefly communicated (Taylor and Wood 1952). 

Other determinations of the neutron total cross sections of a number of 
elements in this same range of energies have been reported at 42 Mev 
(Hildebrand and Leith 1950), at 83 mev (Cook, McMillan, Peterson and 
Sewell 1949), at 95mev (DeJuren and Knable 1950) and at 115 Mev 
(DeJuren and Moyer 1951). Additional measurements of the carbon 
total cross section have been reported (Mott, Guernsey and Nelson 1952, 
Taylor et al. 1951 b). 

The values of the total cross sections presented in this paper are not in 
disagreement with those obtained by other investigators. However, by 
using a narrower energy band than in the work at 42, 83, 95 and 115 Mev, it 
has been possible to look for more detail in the variation of cross section 
with energy. This investigation has shown a dip in the lead Cross section 
at about 60 Mev and has indicated the occurrence of corresponding dips 
for cadmium and copper between 30 and 40 Mev. It has not been possible 
to follow the variation of cross section with energy below 30 Mev since the 
counting rates became very small. . 

Since several parts of the apparatus used in these measurements were 
similar to those described in the previous experiment (Taylor, Pickavance, 
Cassels and Randle 1951 c), the details will not be repeated here. 


* Communicated by the Authors. 
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§2. EXPERIMENTAL 
(i) Principle 
The attenuation of a neutron beam was measured using a ‘ good ’” 
geometry arrangement. The neutron beam was produced by proton 
bombardment of an internal cyclotron target and was then collimated by a 
series of holes in the concrete shielding. There were five sets of holes in the 
two shielding walls (fig. 1, Taylor et al. 1951 ¢), each set corresponding to a 
particular target radius. It was possible to vary the effective energy of the 
neutron beam at a particular radius by decreasing the magnetic field and 
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Variation of the neutron total cross section with energy for lead, cadmium and 
copper in the energy range 14 to 156 Mev. 
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by altermg the target material. For instance the p-n threshold in beryl- 
lium is about 2-0 Mev whereas the threshold in carbon is about 20 Mev. 
A triple coincidence counter telescope detected the neutrons by counting 
the recoil protons emerging from a polyethylene disc placed in the colli- 


mated neutron beam. 


Fig. 2 
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Variation of the neutron total cross section with energy for aluminium and 
carbon in the energy range 14 to 156 Mev. 


The upper energy limit of the neutrons detected was governed by the 
position of the target in the cyclotron, by the value of the magnetic field 
and by the threshold energy required for the production of neutrons by 
protons. A lower energy limit was set by placing absorbers between the 
counters, so that only recoil protons with energy greater than a given value 
were detected. The band of neutron energies used varied from about 
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10 mev at the lower energies of the investigation to about 20 Mev at the 
higher energies. Thus a relatively narrow band of neutron energies was 
used. 
(ii) Geometrical Layout 

The neutron beam, after emerging through a thin window in the 
cyclotron tank, was collimated by a 1 inch diameter hole in a steel bar 
2 feet inlength. This steel bar was inserted in one of five 2 inch diameter 
holes in the main concrete shielding wall, which was 6 feet thick (fig. 1 of 
the earlier report). The 1 inch diameter collimator was placed 26 feet 
from the internal target and the attenuators, 2} inches in diameter, were a 
further 8 feet behind the collimator. Exposure of an x-ray film at the 
position of the attenuators, showed that the beam was closely defined to 
11 inch diameter. The polyethylene disc, 25 inches in diameter, and the 
counter telescope were placed 30 feet behind the attenuators in the shielded 
tunnel. The target, attenuators and polyethylene disc were aligned by 
means of a theodolite. 

(iii) Attenuators 


The attenuators were machined to 24 inches in diameter and the various 
lengths were chosen to correspond to about two mean free paths for the 
particular neutron energy. A remotely controlled sample changer, 
operated from the counting room, held the five attenuators in position and 
also had a position for measurements with the unattenuated beam. The 
sample changer repeatedly positioned the attenuators to an accuracy of 
1/32 inch. Chemically pure materials were used for the attenuators. 


(iv) Neutron Detector 


A triple coincidence counter telescope, using gas counters operated in 
the proportional region, was used to detect recoil protons emerging from 
the polyethylene disc placed in the collimated neutron beam. The 
counters, which were placed at a small angle to the main neutron beam, 
were rectangular in cross section, but otherwise similar to those used in the 
experiment at 153 Mev (‘Taylor et al. 1951 c). A similar arrangement of 
amplifiers, coincidence unit and scalers was used to measure the triple 
coincidences. As previously, random triple coincidences were estimated 
by measuring the coincidences between double coincidence pulses and the 
delayed third channel pulses. The random triple coincidence rate was 
never more than 1°, of the true triple coincidence rate. 

The thickness of the polyethylene disc used at any energy was chosen so 
that it was a small fraction of the total absorber in the telescope. Carbon 
absorbers were placed between the second and third counters in order to 
set a lower energy limit to the recoil protons detected. 


(v) Neutron Monitor 


The intensity of the incident neutron beam was monitored by a single 
boron fluoride counter immersed in a block of paraffin wax and placed 
behind a concrete block in the shielding wall, Tests were made under 
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steady cyclotron conditions of the variation of counting rate with discrim- 


inator voltage and counter voltage, and conditions were set for minimum 
variation. 


§3. Neutron Enercy Banp anp THE ErrecrivE ENERGY 

Only at the highest energies of this investigation was there any appre- 
ciable contribution from the carbon content of the polyethylene, (CH,),,, to 
the detected protons. This arose from using an energy band greater than 
the 13 Mev threshold of the ?C(np)”B reaction. However, even at these 
energies the contribution was less than 3% of that coming from the 
hydrogen content. At the lower energies the contribution was less 
than 1%. 

Neutron Total Cross Sections in 10-27 em? 


Element 


7 3171 +46 
1150 +20 | 1825 -+35 | 2428-442 | 3119 +48 
2234-7-6 | 1100-20} 


3288-150 
892 +22 -|- 50 | 3376 --52 
789 +12 = Lie e 3324+48 
126+3 | 784+5 


571-46 | 1188+11 | 2260-+22| 3341431 
3146461 
508-45 
3083-52 
408+4 | 83348 |1702-+10| 2686437 
390-43 | 799-417] 1627422] 2568+30 


355-43 | 729-411] 1490+20| 2349+19 


330-63. |-677-11)| 1376 +18 


* Previous results with energies corrected. 


Since 97° or more of the protons counted were hydrogen recoils, the 
proton lower energy limit of the telescope corresponded to a definite and 
known neutron energy. By measuring the number of recoil protons as a 
function of absorber thickness and using a range-energy relation, it was 
possible to calculate the effective neutron energy. This measurement also 
provided a check of the maximum neutron energy against the value 
expected from the position of the target in the cyclotron, the value of the 
magnetic field and the threshold energy required for production of neutrons 
by protons. The agreement between the maximum neutron energy 
estimated from the range-energy relation and the expected value gave an 
estimate of the possible error in the absolute value of the effective energy. 
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Thus not only was a narrow energy band of neutrons used, but also the 
effective neutron energy was known within close limits. 

A comparison of the range in carbon and aluminium of the external 
proton beam has suggested that the effective energy of the previous 
experiment quoted as 156-3 Mev may be nearer 153-3 Mev. 


$4, EXPERIMENTAL PROCEDURE 

The experiment was performed by placing each attenuator in turn in 
the neutron beam and comparing the counting rate with that obtained 
when no attenuator was present. This cycle of operations was then 
repeated until the desired statistical accuracy had been obtained. The 
cross sections obtained from each individual cycle showed no greater 
statistical deviation from the weighted mean than that expected from the 
total number of counts recorded. . 

The final result was taken as the weighted mean of the individual cycles. 


§5. RESULTS 

The values of the total cross sections obtained in the neutron energy 
range 30 to 139 Mev for carbon, aluminium, copper, cadmium and lead are 
given in the table together with results previously obtained. 

A small correction has been added to the measured cross section in order 
to take account of those neutrons which were elastically scattered into the 
detector. The lengths of the attenuators and the geometrical arrangement 
were such that this correction was less than the experimental standard 
deviation. 

The cross sections obtained in this investigation, together with those 
obtained by other investigators in the same energy region, are shown in 
the figures. 

§6. Discussion 

The lead cross section shows a dip at intermediate energies, the differ- 
ence between the value at 79 Mev and the miminum value at 60 Mev being 
ten times the standard deviation on either point. Assuming that the 
variation of the lead cross section with energy is correctly depicted by the 
values obtained in this experiment, the lower values shown at 83 and 
95 Mev (Cook et al. 1949, DeJuren and Knable 1950) can be partially 
accounted for by the broader neutron energy band used in those deter- 
minations. 

Including the 25 Mev results (Sherr 1945) and the 14 Mev results 
(Amaldi, Bocciarelli, Cacciapuoti and Trabacchi 1947, Coon 1952) as 
shown in the figures, a similar but slightly smaller dip is indicated in the 
cross section of cadmium and copper between 30 and 40 Mey. Also the 
value of the aluminium cross section at 14 Mev (Coon 1952) is lower than 
the value at 30 Mev by four times the standard deviation on either point, 
again indicating the possibility of a dip below 30 Mev. 

A partial explanation of the shape of these curves is suggested in an 
accompanying letter by Lawson (1953). 
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A Diffraction Effect illustrating the Transparency of N ucler 
to High Energy Neutrons 


By J. D. Lawson 
Atomic Energy Research Establishment, Harwell* 


[Received November 17, 1952] 


Tue transparency of nuclei to high energy neutrons has been strikingly 
demonstrated by recent cross section measurements by Taylor and Wood 
described in a companion paper (1953). 

According to the ‘ continuum’ nuclear model of Bethe (1940), later 
extended by Fernbach, Serber and Taylor (1949), the nucleus can be 
roughly represented by a sphere possessing a complex refractive index 
to the neutron waves. The real part of the refractive index is due to the 
nuclear potential well, and the imaginary part represents a loss of neutrons 
by absorption or inelastic scattering. The value of V, the depth of the 
potential well, is determined by the neutron binding energy and the energy 
of the highest momentum state of the nucleons on the assumption that 
they form a degenerate Fermi gas; it is of order 30 Mev (Fermi 1950). 
The absorption coefficient may be represented by an imaginary potential 
iV’ but its value is not important to the present argument so long as V’ 
is not large compared with V, and is small compared with the incident 
neutron energy L. 

From this model it is easy to find the change of phase of the neutron 
wave if a slab of nuclear matter of depth D is placed with its faces parallel 
to the wavefront. Now if a sphere of diameter D, with a refractive index 
not very different from unity, is introduced instead of a slab there will 
be behind it a circular region of diameter approximately equal to D where 
the average value of the phase retardation is rather less than that which 
would be caused by a slab of thickness D. In the case of the sphere the 
diffraction scattering will be a maximum when this phase change is an 
odd multiple of 7, and a minimum when it is an even multiple. 

The phase change caused by the slab is given by 6=D(k,—k,) where 
k, and ky are the propagation constants inside and outside it. From the 
elementary relation H=}Mv?=h?/2M )? we find ky=27/Ap= 21/27 M 1212 /h 
and ky=2y/27M1'? | H+V-+iV’ |!?/h. Now if V and V’ are small 
compared with H, k,—ky=(21/27M!"H1/h)\(V/2E), so that for the slab 
$/7= (ky —ko) Dl r= 1/2D M1? V hE} = 108D/H1? where D is in cm, and E 
inMev. For the sphere the phase change will be given by $/7= 103%<D/H1? 
where ¢ is a constant rather less than 1. A phase change of nz will 


* Communicated by the Author. 
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therefore occur at an energy given by H=10°%7?D?/n?, Taking rough 
values of 0-6 for « and 2-8 x 10713 41/3 for D, the value of # for which n—1 
is shown in the table for various elements. 


Element C ee ee od) Pb 
Value of # for which d=z 14 26 46 64 100 


It will be seen that these coincide, as nearly as might be expected from 
these very naive considerations, with the peaks in the cross section curves 
as measured by Taylor and Wood. The formula predicts other maxima 
and minima at lower energies than those shown, but at such energies 
the conditions for the validity of these considerations are not expected 
to hold. 


I wish to thank Dr. J. M. Cassels and Mr. A. E. Taylor for illuminating 
discussions on this problem. 
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A New Parameter for Use in Discussions of the Shell Model of the Nucleus 


By N.. Featuer, F.R.S. 
The University, Edinburgh* 


[Received November 4, 1952] 


Iv is well known that the shell model and the liquid-drop model of the 
nucleus derive from approximations to a theory of nuclear structure which 
are distinct and to a large extent mutually contradictory. Each model, 
however, has had considerable success in correlating the results of experi- 
ment in a wide field. It is the object of this note to draw attention to a 
method of using certain experimental results which focuses attention 
on the differences between the properties of real nuclei and those predicted 
by the v. Weizsicker formula (taken as representing the liquid-drop 
model). If the shell model is to be preferred to the other, then obviously 
it must be capable of providing a plausible explanation of these differences. 


A A 
As in a previous note (Feather 1952), we use B, gz): B, Z to 


represent, respectively, the energies of binding of the last proton and 


Pein aes ee 
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A E 
the last neutron in the nucleus ( 7) , and we define a quantity »4,(Z) 


by the relation 
2Z+D 2Z+D—1 
p4,(Z)=8, ( Z, ) By ( 7, ; 


A 
Here D is the isotopic number of the nucleus a and p)4,,(Z) is a measure 
of the extent to which the energy of binding of the last (Z*") proton in this 
. (A—1] : 
nucleus exceeds that of the last proton in ( Z ) . Using a result 


established in the above-mentioned note (Feather 1952), we also have 
2Z+D 2Z+D—1\ _ 2Z+D\ _ 2Z+D—1 
B, ( Z, ) —B, ( Z =B,, Z 8 Treat ’ 
so that p4,,(Z) likewise measures the extent to which the energy of binding 


of the last ((Z-+D)*®) neutron in the nucleus i exceeds that of the 


(Z+D)*® neutron in The parameter of interest for the discussion 


A—l 
(7-1) 
of nuclear models is now »4’,, (Z), where 

p4’,(Z)=p-24,(4+1)—p4,(Z). 
It will be observed that p4’,,(Z) represents the increase, for unit increase in 
the proton number Z (from Z to Z+1), in the difference between the 
energies of binding of the last proton in nuclei of mass numbers 4 and A —1, 
respectively (as also the corresponding quantity in respect of neutron 
binding energies). Other formal properties of p4’,(Z) are as follows : 


(i) If 7,- a represents the total energy of negative £-disintegration, 


ground-state to ground-state, for the species , then 


A 
Z 
pd ,(Z)=Py- (aera) a7, (irate ay pers) 


Z Z—1 
_ (224-D--2 
ay (EP), 
Here we may note that the four species involved in this relation comprise 
one example each of the neutron-proton parity types e—e, 0—0, €—0, o—e, 
whatever the parities of Z and D. 

(ii) Assuming the four-constant v. Weizsaicker formula for the ‘ exact ’ 
nuclear mass (that is omitting any spin or parity-type term) the terms in 
a and y (‘ volume’ and ‘surface ’ effects) vanish identically in the 
general expression for p4’,(Z)—and the term in B (‘ asymmetry ’ effect) 
vanishes identically when D—1. 

(iii) On the same assumption, and with accepted values for the vy. 
Weizsicker constants «, B, y and 8, for all values of D the numerical 
magnitude of p4’,(Z) is negligibly small for all but the very smallest 
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values of Z (e.g. for D=2, Z=7; D=3, Z=16; D=51, Z=90 calculated 
values of »A’,,(Z) are, respectively, 0-052, 0-0096 and 0-0008 Mev). 

Since p)4’,(Z) can rarely be known from experiment to better than 
+0-05 Mev, it will be seen from (iii) that experimental values of this 
parameter significantly different from zero point to a failure of the 
v. Weizsacker liquid-drop model formula completely to account for the 
facts in individual cases—and it will be seen that this conclusion is 
effectively valid over the whole range of D and Z. 

A survey of experimental data leading to the calculation of more than 
600 neutron and proton binding energies has recently been completed, and 
an account will shortly be published with full discussion using the para- 
meter »4’,(Z) and other criteria. For the purposes of this note it is 
sufficient to state that the ‘experimental’ curves of p4’,(Z) against 
Z (D constant) show interesting, and obviously significant, features. 
Thus the 5 points which are available for D=1 can be represented 
within the limits of experimental uncertainty by the straight line 
A’=0, but the 12 points available for D=2 represent positive values 
of A’ ranging from 0-4 to 4:0 Mev. Experimental values of J’ consider- 
ably different from zero are less numerous as D (and Z) increases—thus 
the 9 points for D=47 fall on the line 4’=0 almost within the limits of 
error—but variations in 4’ (D constant) as large as 1-5 Mev certainly 
persist as far as D=15. 
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The Surface of Beryl and Cadmium Iodide Crystals 


By H. F. Kay and BarBaRa APPELBE 
H. H. Wills Physical Laboratory, University of Bristol* 


[Received November 4, 1952] 


OPTICAL microscope studies of the prism faces of beryl crystals (Griffin 
1951) revealed unimolecular 7:9 A steps of surprisingly high visibility 
which nevertheless varied considerably over different areas of the crystal 
surface. This strongly suggested the presence of surface features peculiar 
to the step, which, were just smaller than the resolution of the optical 
apparatus employed. The electron microscope was employed to investi- 
gate these surface features by standard replica techniques. 
The prism faces (1010) of beryl crystals were found to be covered with 
small pits which are attributed to some etching process having occurred. 
subsequent to growth. The high visibility of the steps is ous to the 
higher concentration of pits along the step edges. Figure 1 (figs. 1, 2, 


RG ae 
* Communicated by the Authors. 


106 Correspondence 


Plate 1, figs. 3, 4, Plate 2) is an electron micrograph of a single-stage 
shadowed formvar plastic replica taken from a prism face, and shows this 
concentration of the pits along lines representing the step-edges. The 
variation in optical visibility results from changes in size (2 u to 0-05 1), 
depth (0-4 1 to 0-02 u) and distribution of the pits. ne 

In general, no characteristic regularity of shape or linking together 
seems to occur, for quite different types of corrosion can be seen even on 
different areas of the same face, but higher resolution two-stage methyl- 
methacrylate silica replicas of a frequently occurring type of pit show 
that here at least they appear to be of more definite shape and are striated 
parallel to the basal plane at intervals of ~ 350-400 A. Figures 2 and 
3 show two examples of this kind. 

Replicas of the terminal faces (0001) of beryl revealed multimolecular 
steps and the absence of the general pitting such as occurred on the 
prism faces. Mottled patches approximately a few microns in extent 
were found to be scattered over the terminal surface and these are 
possibly finely pitted areas, the pits themselves being on the limit of the 
horizonal resolution of the plastic replicas. The absence of pits along the 
step edges on the terminal faces explains the paradox, observed optically 
(Griffin 1951), that the visibility of the multi-molecular steps is much 
less than that of the unimolecular steps on the prism faces. 

The growth of cadmium iodide crystals as revealed by the optical 
microscopy (Forty 1952) is consistent with the spiral mechanism predicted 
by the dislocation theory. Electron optical observations were under- 
taken on these crystals to try and elucidate the nature of the large 
dislocation sources observed. Owing to their high molecular weight, the 
crystals could not be examined directly by electron transmission, and 
replica techniques were employed to study the surface features. The 
small cadmium iodide crystals were grown between microscope slide and 
a cover-slip from a drop of warm saturated solution, A polystyrene in 
benzene replica was then taken of the crystals on the slide. Although 
polystyrene replicas were found to be the strongest, they almost invariably 
broke along the step edges due to the fact that the step heights, usually 
up to 2000 4, were much greater than the thickness of the replica. 
Figure 4 shows a crystal whose growth is activated by two dislocations 
of opposite hand and of multiple strength, the centres of which are 
marked by hexagonally shaped prismatic holes approximately 44 in 
diameter, and having three of their edges parallel to the growth steps. 
(See inset.) Forty (1952) has seen indications of such hollow dislocation 
centres in optical photographs of cadmium iodide and they have been 
predicted on purely theoretical grounds by Frank (1951) 


. 
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Fig. 1 
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The surface of a prism face of a beryl crystal. 
with gold-palladium at 10°. 


shadowed with gold-palladium at 10°. 


seen on the prism face of beryi crystals. 
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108 Notices of New Books and Periodicals received 


Econometrics. By G. TrytNER. (Chapman and Hall.) Price 46s. 


Tuts book will be of value to those who know both a little economics and a 
little mathematical statistics, although the standard of exposition is not as 
high as the subject deserves. The first part of the book is introductory. 
The second deals with multivariate analysis. The last, modestly entitled 
‘ Some topics in time series analysis ’, efficiently treats problems of importance 
to economists which are not discussed in other text-books, and is therefore 
of great value to the student or the lecturer. ge 

Throughout the book the examples are apposite and helpful, while the 
footnotes would make a full and useful bibliography—the absence of which 
is to be regretted. He AST 


A Short Table of the Bessel Functions In+4(#), (2/7)Kn+3(xv). Prepared by 
C. W. JonEs on behalf of the Mathematical Tables Committee of the Royal 
Society. (Cambridge University Press.) Price 6s. 6d. 

Tae functions tabulated here are, with the corresponding ranges of argument 

and order, 


a-"—4],, 4 (2) v=0(0-1)5 ; n=0(1)10 
(2/r)a"+*# Kn 4. 4(2) a=0(0-1)5 ; n=-0(1)10 
e- In 4 4(a) a=5(0-1)10 ; n=0(1)10 
(2/r)e*Kn 44(2) v=5(0-1)10 ; n=0(1)10 


which for obvious reasons are more convenient to tabulate than the Bessel 
functions themselves. Seven or eight significant figures are given almost 
everywhere and modified second differences are provided 


Metallurgical Equilibrium Diagrams. By W. Hume—-Rotuery, J. W. 
CuristTiAn and W. B. Pearson. (The Institute of Physics, London.) [pp. 310.] 
Price 50s. 


Tals is a good book ; but the title may be a little misleading. It deals primarily 
with the experimental methods employed on establishing the diagrams and 
not at all with the results of such investigations. There is in addition a 
brief theoretical introduction. The authors have interpreted these terms of 
reference generously and have ranged far in search of material. There are large 
sections on thermometry and thermal analysis and less extensive chapters on 
x-ray, dilatometric and magnetic methods, and on the techniques of metallo- 
graphy. The book includes 239 diagrams, 202 references to the literature and 
a copious index. Although intended primarily for research students in 
metallurgy, it will prove valuable to all metallurgists and metal physicists and 
indeed to all workers concerned with laboratory experiments at high tempera- 
tures. 


Progress in Metal Physics, Vol. 111. By B. Cuaumers (Editor). (London ; 
Pergamon Press Ltd., 1952.) [Pp. 334.] Price 48s. 


Tae high standard of previous volumes is maintained in this latest group of 
reviews of progress in certain fields of metal physics. The subjects covered 
in this volume are crystallography of transformations, low temperature 
properties, electron theory, twinning, ferromagnetism, x-ray strain measure- 
ments, re-crystallization, and grain boundaries. Ae 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


